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Measurements  of  electromotive  force  (EMF)  and  the  iso- 
piestic  technique  were  utilized  in  the  determination  of 
factors  which  influence  the  thermodynamic  activity  and 
related  properties  of  individual  species  in  complex  saline 
mixtures. 

EMF  measurements  in  cells  without  liquid  junction  were 
carried  out  to  determine  the  mean  activity  coefficient  of 
HC1,   y+>  in  electrolyte  mixtures  resembling  seawater.   In 
mixtures  of  HCl-MgCl2,  it  was  found  that  y+  decreases  in  a 
linear  fashion  as  HC1  is  replaced  by  MgCl„  while  maintaining 
constant  total  ionic  strength  and,  in  this  respect,  obeys 


IX 


Harned's  Rule.   These  measurements  were  carried  out  at  four 
ionic  strengths,  namely  I  =  0.1000,  0.3809,  0.6729,  and 
0.8720  mol  kg-1  at  nine  temperatures  from  5  to  45°C.   The 
three  higher  ionic  strengths  correspond  to  seawater  of 
salinities  20,  35,  and  45°/00.   In  addition,  studies  were 
made  on  the  system  HCl-NaCl-MgCl2  at  I  •=  0.6729,  the  molal 
ratio  of  NaCl  to  MgCl2  being  maintained  at  7.202  as  in 
natural  seawater.   The  Harned  coefficients  obtained  from 

all  measurements  were  found  to  decrease  slowly  with  tempera- 

tr 
ture.   The  trace  activity  coefficients  of  HC1,  y±  ,  were 

calculated  from  the  Harned  equation  in  two  ways.   First,  by 
combining  Harned  coefficients  from  the  HCl-MgCl2  study  with 
those  from  an  earlier  HCl-NaCl  study,  weighting  the  two  with 
respect  to  their  relative  ionic  strength  contributions  to 
seawater.   Secondly,  they  were  calculated  using  Harned  coef- 
ficients from  the  HCl-NaCl-MgCl2  study.   All  values  were 
found  to  be  identical  with  those  measured  earlier  in  total 
artificial  seawater  (excluding  sulfate)  at  all  ionic  strengths 
and  temperatures  included  in  the  experiments. 

The  concept  of  single-ion  activity  and  the  problem  of 
liquid-junction  potentials  in  pH  measurement  were  addressed 
in  a  study  of  pH  in  synthetic  seawater  and  saline  media  simu- 
lating serum.   Measurements  of  EMF  in  a  cell  with  liquid 
junction  at  25°C  were  carried  out  in  buffered  artificial  sea- 
water at  salinities  of  30,  35,  and  40°/00,  each  containing  a 
Tris  buffer  at  molalities  of  0.02,  0.04,  and  0.06. 


Comparison  of  the  pnu  values  calculated  by  the  operational 
definition  with  those  determined  from  measurements  in  cells 
without  liquid  junction  demonstrates  that  a  residual  liquid- 
junction  potential  exists  which  varies  with  salinity.   Within 
a  given  salinity,  however,  pm„  values  calculated  by  the 
operational  definition  are  identical  regardless  of  Tris 
buffer  concentration. 

Evidence  is  presented  that  pH  measurement  of  blood  plasma 
and  other  clinical  media  at  ionic  strengths  of  0.16m  vs.  the 
NBS  standards  may  involve  residual  liquid-junction  errors 
amounting  to  0.03  to  0.05  pH  unit.   EMF  measurements  of  cells 
both  with  and  without  liquid  junction  indicated  that  residual 
liquid-junction  effects  may  be  nearly  eliminated  by  matching 
the  ionic  strength  of  the  standard  to  that  of  the  sample.   A 
phosphate  buffer  composed  of  0.005217m  KH2P04,  0.018258m 
Na„HP0   and  0.1m  NaCl  (buffer  ratio  1:3.5  and  I  =  0.16m)  is 
to  be  preferred  to  the  NBS  phosphate  buffer  (I  =  0.1m)  for 
measurements  of  this  sort. 

Isopiestic  measurements  were  carried  out  at  25°C  on  mix- 
tures of  electrolytes  having  to  do  with  seawater,  NaCl-SrCl2 
and  NaCl-Na2C03,  and  the  resulting  osmotic  coefficients  were 
determined  in  the  ionic  strength  range  0.3  to  6.0m.   The  data 
were  fitted  to  the  current  treatments  of  Pitzer  and  Scatchard 
for  solutions  of  mixed  electrolytes.   Activity  coefficients 
were  calculated  for  the  individual  solutes  in  each  mixture 
using  the  best-fit  ion  interaction  parameters  from  the  Pitzer 


treatment.   These  are  discussed  in  light  of  the  excess  free 
energies  of  mixing  calculated  from  the  corresponding 
Scatchard  parameters. 
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CHAPTER  1 
INTRODUCTION 


When  a  salt  such,  as  sodium  chloride  is  dissolved  in 
water,  the  resulting  solution  exhibits  a  comparatively  high 
electrical  conductivity  showing  that  charged  ions  are  pres- 
ent.   The  thermodynamic  properties  of  the  solution  may,  in 
part,  be  described  by  electrostatic  interactions  of  the  ions 
with  each  other  and  by  interactions  of  ions  with  water  mole- 
cules.  The  concentration  of  the  dissolved  species  in  the 
solvent,  and  thus  the  magnitude  of  these  effects,  will  have 
a  major  influence  on  the  thermodynamic  activity  and  related 

properties. 

1  —  7 
Classical  thermodynamics    gives  a  general  formula  for 

the  chemical  potential  of  any  solute,  i,  in  an  ideal  solution 
as : 

ui  =  u|  +  rt  In  x±  Cl.l) 

where  x^^  is  the  mole  fraction  of  the  solute  and  y°  is  its 
chemical  potential  in  the  standard  state,  i.e.,  when  x.  is 
unity.   The  temperature  is  given  by  T,  and  R  is  the  gas 
constant.   In  an  electrolyte  solution,  however,  coulombic 
ion-ion  forces  are  present  and  any  attempt  to  measure  the 
chemical  potential  would  result  in  a  departure  from  equation 


1.1,  the  magnitude  of  which  is  a  measure  of  the  aoa- ideality- 
associated  with  these  interactions. 

In  order  to  correct  equation  1.1,  an  empirical  correc- 
tion factor  may  be  applied  to  modify  the  concentration  term: 

(i.  =  U'  +  RT  lax.f.  (1.2) 

This  effective  concentration  term,  x.f^  ^  a^  is  termed  the 
activity  of  the  species  i,  and  f .  is  the  "rational"  activity- 
coefficient.   The  concentration,  may  be  just  as  easily  ex- 
pressed in  molality,  m. ,  where  the  activity  coefficient  must 
take  on  a  different  value  normally  designated  by  the  symbol 
y..   From  equation  1.2,  it  can  be  seen  that  the  activity  must 
be  unity  in  the  standard  state.   Thus,  the  standard  state  on 
the  molality  scale  is  that  of  a  hypothetical  solution  where 
the  mean  ionic  molality  and  activity  coefficient  are  unity, 
but  with  the  reference  state  chosen  such  that  a./m^  =  yi  =  1 
when  m.  approaches  zero. 

The  change  in  the  Gibbs  free  energy  of  any  thermodynamic 
system  as  a  function  of  temperature,  pressure,  and  composi- 
tion is  given  as: 

dG  =  -SdT  +  VdP  +  £uidni  11.3) 

i 

Therefore,  mathematically,  the  chemical  potential  of  a  species 

i  is  the  partial  derivative  of  the  free  energy  with  respect 

to  moles  of  i  when  temperature,  pressure,  and  moles  of  other 

substances  remain  constant;  it  is  often  written  as  G.. 

5i  =fe  J  CI. 4) 

1   \5ni/T,P,n. 
3 


To  measure  this  quantity  in  an  ionic  solution,  one  would 
need  to  measure  the  change  in  free  energy  with  change  in 
concentration  of  one  ionic  species  only.   Obviously,  this 
would  be  empirically  difficult  if  not  impossible.   One  can 
only,  therefore,  measure  the  activity  coefficient  of  a  net 
electrolyte,  that  is  at  least  two  ionic  species  resulting 
from  a  single  neutral  electrolyte. 

If  one  considers  a  symmetrical  electrolyte  of  the  1:1 
charge  type,  such  as  EC1 ,  the  sum  of  the  chemical  potentials 
of  the  resulting  ions  as  a  function  of  molality  may  be  given 
as : 

GH+  +  GC1_  =  (G°+  +  G^.)  +  vRT  In  (m^+rn^-)1^ 

+  vRT  In  (Yh+^c1-)1/V       (1'5) 
where  v  =  v   +  v~  is  the  number  of  ions  per  molecule  of  elec- 
trolyte and  in  this  case  is  equal  to  2.   Defining  the  respec- 
tive mean  values  as: 

GH+  +  GC1- 

b^  =  - 


G?  = 


2 


G°+  +  G°._ 


,   v+  v-  Nl/v 
m±  =  ^raH+mCl-^ 

,    v+  v-  sl/v 
Y±  =  (YH+YC1-) 

an  expression  can  be  written  for  the  mean  chemical  potential 
which  is,  in  effect,  the  average  contribution  of  a  mole  of 
ions  to  the  free  energy  of  the  system. 

G+  =  G°  +  vRT  In  m+  +  vRT  In  y+  CI -6) 


One  now  has  a  relation  between  the  free  energy  of  a  state  in 
which,  the  ions  are  infinitely  far  apart  and  a  state  corres- 
ponding to  a  given  concentration  and  its  mean  activity  coeffi- 
cient, y+j  which  is  experimentally  determinable..   To  split 
this  value  of  y  into  individual  ionic  activity  coefficients 
will  require  a  theoretical  and  partly  arbitrary  approach  to 
be  discussed  later. 

Another  way  in  which,  the  non-ideal  behavior  of  an  elec- 
trolyte solution  may  be  described  is  by  means  of  the  osmotic 
coefficient ,  which  may  be  regarded  as  a  measure  of  the  nonr 
ideality  of  the  solvent  due  to  the  presence  of  ions.   It  is 
defined  in  such  a  way  as  to  become  unity  in  the  pure  solvent. 
A  rational  activity  coefficient  for  the  solvent  would  fail  to 
emphasize  this  adequately  since  the  mole  fraction  of  the  sol- 
vent is  normally  so  large,  f  would  never  deviate  much-  from 
unity.   For  this  reason,  the  "rational"  osmotic  coefficient, 
g,  is  defined  by: 

/    vmW  \ 
in  as  =  g  in  xs  =  -g  In  ^1  +  ^ )  U-7). 

where  W  is  the  molecular  weight  of  the  solvent ,  a  is  its 
s  s 

activity,  and  m  is  molality.   The  practical  osmotic  coeffi- 
cient, 4> ,  on  the  molality  scale  is  defined  by: 

-vmW 
ln  as  =  Toof  *  Cl. 8) 

Therefore,  for  an  aqueous  solution  for  a  single  electrolyte, 

equation  1.8  becomes: 

-55.51    In   aw  =   vm<j>  Cl<9) 


In  very  dilute  solutions,  the  mole  fraction  of  the  solute, 
x.,  may  be  described  by: 

n .     n. 

xi =  h— —  B  IT  <1.10) 

w    1    w 
where  n^  and  nw  are  the  number  of  moles  of  solute  and  water, 
respectively.   Therefore  from  equation  1.9  one  finds: 

-m  In  a    -In  a 

«*  =  ^-J!L~-    "i^  Ci.U) 

When  x.  is  very  small,  the  last  term  may  be  approximated  by: 


-In   a 
w 

= 

-In  a  In  a 
w                 w 

x  . 

1 

1  -  x  In  x 
w                 w 

(1.12) 

From  this  it  can  be  seen  that  the  osmotic  coefficient  is  approx- 
imately proportional  to  the  ratio  of  the  logarithms  of  activity 
to  mole  fraction  and,  in  this  respect,  behaves  much  like  the 
activity  coefficient  where  f.  =  a./x.. 

The  relationship  between  the  rational  and  the  molal  activ- 

4 
lty  coefficient  can  be  shown  to  be: 

f±   =   Y+Cl  +  0.001vWgm)  (1.13) 

and  that  between  the  rational  and  molal  osmotic  coefficient  as: 

vmW 
g  ln  xs  =  Toot  *  Cl. 14) 

In  the  case  of  any  system  at  equilibrium  and  maintained 

at  constant  temperature  and  pressure,  the  sum  of  the  chemical 

potentials  of  all  species  must  be  zero;  as  shown  by  the  Gibbs- 

Duhem  relation, 

VSA  +  VS  +  V^C  +  '••  =  °  CI.  15) 


where  n.  denotes  the  number  of  moles  of  the  ith  species. 
Therefore,  for  a  solution  containing  only  one  solute  in 
water: 

n.dG.  =  -n  dG  (1.16) 

11     ww  v -     ' 

The  Gibbs-Duhem  relation    is  important    in   relating   the   osmotic 
and  activity  coefficients.      Recalling  that: 

G.    =   G°    +   vRT   In   a. 

11  l 

and   G     =   G°    +   RT   In   a 

WW  w 

one   can   express  the  Gibbs-Duhem  relation   as: 

vm   din    Cr+ml   =   -55.51   din    a  CL17) 

Combining  this   result   with  equation  1.9,    the  osmotic   and 
activity   coefficients  may  be  related  by: 

C*    -   l)^jp  +   d<J)   =   din  Y±  (1.18) 

and  upon  integration  this  may  be  shown  to  give: 

m 
4=1+  1/m  /  m  din  y+  (1.19) 

0 

Alternatively,  activity  coefficients  can  be  calculated 
from  osmotic  coefficients  by  expressing  the  above  relation- 
ship in  the  form: 

m 
In   y+   =    C<j>    -   1)    +   /    C<J>    -   1)    din  m  (1.20) 

0 


CHAPTER  2 
ACTIVITY  COEFFICIENTS 


Theoretical  Background 

Electromotive  force  measurements  have  been,  useful  in 
providing  a  large  part  of  our  existing  t he rmo dynamic  data 
for  electrolyte  solutions.  They  serve  as  the  most  common 
and  accurate  means  by  which  to  view  ionic  activity. 

The  reaction  in  a  galvanic  cell  is  one  in  which,  oxida- 
tion and  reduction  occur  simultaneously;  that  is,  electrons 
are  removed  from  one  species  wh-ich.  is  said  to  be  oxidized 
and  supplied  to  another  species,  said  to  be  reduced.   The 
reaction  may  be  regarded  as  the  sum  of  these  two  half-reactions, 
each  of  which  occurs  at  a  respective  electrode.   If  the  two 
electrodes  are  connected  to  a  voltmeter  with  a  high,  input 
resistance,  such  that  no  reaction  is  allowed  to  occur,  the 
resulting  voltage,  E,  is  a  measure  of  the  potential  for  reac- 
tion.  This  potential  is  related  to  the  molar  free  energy 
change,  AG,  for  the  reaction  by: 

AG  =  -nFE  (2.1) 

where  F  is  the  Faraday  (.96,487  coulombs/mole).,  n  is  the  num- 
ber of  moles  of  electrons  exchanged,  and  E  is  in  volts.  The 
Gibbs  energy  change  is  then  in  joules  and  is  a  measure  of  the 


maximum  work  that  the  cell  reaction  can  produce,   A  nega- 
tive AG  indicates  a  spontaneous  reaction,  while  a  positive 
AG  indicates  the  reaction  will  not  occur  spontaneously. 

Considering  the  following  general  chemical  reaction: 

aA  +  bB  +  . , .  =  uU  +  vV  +  . . .  (2 . 2 ) 

then  the  change  in  free  energy  which  accompanies  the  reaction 
is  the  sum  of  all  the  chemical  potentials  of  the  products 
times  their  respective  number  of  moles  minus  the  same  quan- 
tities for  the  reactants: 

AG  =  uGn  +  vGy  +  . , ,  -  aGA  -  bG   -  . , ,         C2 . 3 ) 

If  one  writes  for  AG0 : 

AG0  =  uG°  +  vG°  +  ...  -  aG°  -  bG°  -  ,,,        C2.4} 

and  recalls  equation  1,1,  then: 

u  v 
a^ay  .. .  . 

AG  =  AG0  +  RT  In  -^-4 C2.5). 

aAaB  * • • 

which  is  the  general  form  for  the  free  energy  of  a  reaction. 
Combining  relation  2.1  with  the  above  gives  the  Nernst  equa- 
tion : 

u  v 
E  =  E.o  _  R|     *U*V  C2.6) 

nE     a  b  u  '  -* 

aAaB  ' * ' 

which  relates  EMF  to  the  activities  of  the  various  ions. 

This  equation  will  be  used  extensively  in  this  work.   The 

term  E°  is  the  standard  potential;  it  is  equal  to  E  when  all 

ions  are  in  their  standard  states. 

Considering  the  following  galvanic  cell  *    containing 

aqueous  HC1 : 
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Pt;  H2Cg,  latm)|HClCm)|AgCl;  Ag(s)  C2.7) 

the  cell  reaction  may  be  regarded  as  a  combination  of  two 
half  reactions: 

fH2  =  H+  +  e"  C2.8) 

and  AgCl  +  e~  =  Ag  +  Cl~  (,2.9) 

which,  when  combined,  give  the  total  cell  reaction: 

4H2(g,  latm)  +  AgCl(s)  =  AgCsO  +  H?m)  +  Cl(m)      C2.10) 

The  cell  2.7  is  written  such  that  oxidation  is  taking  place 
at  the  left  electrode,  the  anode,  and  reduction  takes  place 
at  the  right,  the  cathode.   A  positive  EMF,  or  negative  AG, 
indicates  that  the  reaction  will  occur  spontaneously  as 
written.   The  EMF  of  the  cell  is  given  as: 

E  =  E°  -  ~  In  aH+acl^  C2.ll) 

since  hydrogen,  silver,  and  silver  chloride  are  all  in  their 
standard  states.  This  may  be  expanded  to  give  E  in  terms  of 
the  molalities  and  mean  activity  coefficients  of  the  respec- 
tive species  as  follows: 

E  =  E°  -  Sj£  In  mH+mcl_  -  ^  In  y2±  C2.12) 

Therefore,  it  can  be  seen  that  EMF  measurements  are  useful  in 
that  they  easily  yield  the  mean  activity  coefficient  or  log  yH 

logy±  =  -H(E  -  E°)/k  +  logmH+mcl_J         C2.13) 
where  k  =  2.303RT/nF. 
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Solutions  of  Single  Electrolytes 

A  theoretical  treatment  of  the  thermodynamic  properties 
of  electrolyte  solutions  must  deal  with  both  long-range  inter- 
ionic  forces  and  short-range  ion-solvent  interactions.   The 
former  net  effect  of  attractions  and  repulsions  will  tend  to 
decrease  the  free  energy  of  any  particular  charged  species, 
manifested  in  a  decrease  in  the  activity  coefficient.   Also 
stabilization  of  ions  by  solvent  molecules  tend  to  lower 
their  free  energies.   On  the  other  hand,  at  higher  concentra- 
tion ion-solvent  interactions  tend  to  lower  the  vapor  pressure 
of  the  solvent,  and  thus  its  activity,  appearing  to  increase 
the  activity  coefficient  of  the  solute  species.   There  is 
normally  a  concentration  where  these  two  phenomena  are  in 
balance,  resulting  in  an  activity  coefficient  minimum.   How- 
ever, in  dilute  solutions,  interionic  forces  may  be  assumed 
to  dominate. 

The  contribution  of  electrical  interactions  with  other 
ions  to  the  free  energy  of  a  specified  mole  of  j-ions  is 
shown   to  be: 

2  2,T 

_        -z .  e  N     k 

AG  Cel)  =  -J (2.14) 

J  2e    1  +  <a 

This  arises  from  the  well-known  Debye-Huckel  formula  for  the 

time  average  potential  for  an  ion  of  diameter,  a,  with  a 

charge  z  e,  in  the  absence  of  other  external  forces.   The 
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term,  e,  is  the  dielectric  constant  of  the  medium  in  which 
the  ions  are  immersed,  N  is  Avogadro's  number;  k  is  defined 
by 

(i^4)1/2^ 

and  is  obviously  a  function  of  concentration,  ionic  charge, 
temperature,  and  dielectric  constant.   It  is  often  called 
the  Debye-Hiickel  reciprocal  length  since  it  is  this  distance, 
k~  ,  at  which  all  of  the  charge  of  the  ionic  atmosphere  sur- 
rounding a  central  ion  may  be  considered  to  be  concentrated. 
The  Boltzmann  constant,  k,  arises  from  the  assumed  distribu- 
tion function. 

Assuming  that  an  ionic  solution  would  behave  ideally  in 

the  absence  of  these  interionic  forces,  then  the  partial  free 

4 
energy  of  a  mole  of  j-ions  may  be  split  into  : 

G.  =  G. (ideal)  +  AG. (el)  (2.16) 

3  3  3 

The  term  AG. (el)  is  essentially  the  non-ideal  contribution 

containing  the  activity  coefficient,  hence: 

—  2  2 

AG. (el)    -z.e     < 

In  f  .  =  — ^ =   -J (2.17) 

3  RT       2ekT  1  +  ica 

For  the  mean  activity  coefficient,  f + ,  which  is  experimentally 

determinable : 

2  2^,2 

-e     k    v-z-  +  v  z 

In  f+  = -±-± £-£  (2.18) 

2ekT  1  +  Ka   vi  +  vo 
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The  condition  of  electroneutrality ,  v-z-  =  -v„z9,  allows  the 
following  transformation: 

-|z  z  |e     k 
In  f+  =  -±-A (2.19) 

2ekT    1  +  <a 
Substituting  for  k,  equation  2.19  can  be  written  in  decadic 
logarithmic  form  as: 

logf   =   i_f (2.20) 

1  +  Ba/T 

where  A  and  B  are  constants  dependent  on  T  and  e.   The  ion- 

o 

size,  a,  is  given  in  angstroms  for  numerical  convenience. 
A  .  1.8246  x  106 

B  x  108  -  5-029  ^/2109  (2.22, 

At  very  low  concentrations  of  electrolytes,  that  is  at  low  /T, 
the  term  Ba/T  becomes  negligible  compared  to  unity  and  equa- 
tion 2.20  will  take  the  form: 

logf±  =  -A|z1z2|/I  (2.23) 

which  is  known  as  the  Debye-Hiickel  limiting  law.   This  form 
is  useful  in  predicting  activity  coefficient  behavior  in  very 
dilute  solutions.   Therefore  one  could  say,  with  reasonable 
justification,  that  the  numerator  in  equation  2.20  is  the 
effect  due  to  long-range  coulombic  forces,  while  the  denomi- 
nator brings  in  modifications  for  short-range  interactions. 
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By  choosing  3.04  A  as  a  reasonable  value  of  a  for  all  elec- 

12 
trolytes  at  25°C,  Giintelberg    expressed  equation  2.20 

simply  as: 

logf+  =  i_i C2.24) 

1  +  /T 

which,  is  superior  to  the  limiting  law,  but  not  adequate  above 

13 
I  =  0.1m.   Guggenheim,    by  adding  an  adjustable  parameter 

linear  in  concentration,  improved  the  equation  still  further: 

-A|z-zJ/T 

logf+  =  — +  bl  C2.25) 

1  +  /T 


Solutions  of  Mixed  Electrolytes 

Studies  of  mixed  electrolytes  are  important   ~   from 
the  standpoint  that  they  provide  an  experimental  criterion 
for  judging  theoretical  predictions.   From  the  bodies  of 
water  in  the  oceans  to  the  biological  fluids  in  the  body, 

solutions  of  mixed  electrolytes  are  of  great  importance. 

19 
In  1922,  Br^nsted   proposed  his  theory  of  specific 

interaction  which  states  that,  "in  a  dilute  solution  of  con- 
stant total  concentration,  ions  will  be  uniformly  influenced 
by  ions  of  their  own  sign  and  specific  effects  are  to  be 
sought  in  interactions  between  oppositely  charged  ions." 
Employing  a  modification  of  the  Debye-Htickel  expression 
(equation  2.25)  with  a  term,  b,  linear  in  concentration, 
Guggenheim  continued  to  build  a  theory  of  mixed  electrolyte 
solutions.   He  introduced  specific  interaction  coefficients, 
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b,  into  equations  for  the  activity  coefficient  of  an  elec- 
trolyte, B,  in  the  presence  of  another  electrolyte,  C: 

In  yb  =  "1,173/T  +  [2xbMX  +  Cb^  +  bMY)(l  -  x)]m    (2.26) 
1  +  /I 

+      -  +      - 

Here,  M  and  X  are  the  ions  from  B  and  N  and  Y  are  the 

ions  of  C,  both  being  1:1  electrolytes  in  aqueous  solution 

at  25°C.   The  total  molality  is  m,  while  xm  and  CI  -  x)m  are 

the  molalities  of  B  and  C  respectively.   This  form,  taken 

4 
from  Robinson  and  Stokes,   describes  the  activity  coefficient 

and  concentration  on  the  molal  scale.   A  similar  equation  may 

be  written  for  the  electrolyte  C: 

In   yc  =   -1-173vl  +    [2 CI   -  x)bNY  +    (b^  +   bMY)x]m        C2.27) 
1   +    /I 

Now,  letting  y'  refer  to  the  term  in  brackets  on  the  right  of 
equations  2.28  and  2.27,  then  when  x  =  0  and  only  component  C 
is   in   solution   one   can   define: 

ln   YCO)B   =    CbMY   +   bNX)m  ^2.28) 

In   YcC0)    =   2bNym  C2.29) 
In   contrast,    when  x  =   1: 

ln   YB(0)    =    2bMXm  C2'30) 

ln   YCO)C   =    CbNX   +   bMY)m  ^2.31) 
Now  it    can   be  seen   that : 

ln   yB   -    In   Y(0)B   +    Un   YB(0)    -    In   YC0)B)x  (2.32) 

and  similarly  for  C: 

ln  yc  =  In  Y(0)c  +  Cln  Yc(0)  -  In  Y(0)C>C1  -  x)    (2.33) 
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Thus,  the  logarithm  of  the  activity  coefficient  of  either 
component  in  a  mixture  maintained  at  constant  total  molality 

varies  linearly  with  composition. 

20—23 
Earned  and  co-workers      have  made  numerous  EMF 

measurements  in  solutions  of  two  components  at  constant 

ionic  strength,  and  from  their  results  has  emerged  Earned' s 

Rule: 

logTB  =  logy°  -  aBmc  (2.34) 

which  states  that  the  change  in,  logy'g  from  its  value  in  a  pure 
solution  of  B  is  linear  in  molality  of  the  other  component. 
The  Earned  slope,  a„,  is  itself  a  function  of  ionic  strength. 
This  rule  has  also  been  found  to  be  valid  for  mixtures  of 
non-symmetrical  electrolytes  as  long  as  the  total  ionic 
strength  is  maintained  constant.   In  this  case,  the  decrease 
in  logy  of  one  component  is  a  function  of  the  ionic  strength 
fraction,  y,  of  the  other  component. 

EC1  in  Seawater 

In  a  mixture  such  as  seawater,  composed  of  many  dissolved 
electrolytes,  one  could  imagine  the  decrease  in  logy  of  any 
single  electrolyte  to  be  a  function  of  the  respective  contri- 


butions of  the  other  electrolytes  to  the  total  ionic  strength, 

24 
Macaskill  et  al.    have  made  measurements  of  the  cell: 

Pt;  H2(g,  latm)|HCl(mA),  NaCl(m^) | AgCl ;  Ag       (2.35) 
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where  the  molalities  of  the  two  components  were  varied  while 
maintaining  a  constant  total  ionic  strength.   The  activity 
coeiiicient  of  HC1  can  be  calculated  from: 

logyHC1  =-HE  -  B-/k  +  logmHC1(.mHC1  +  m^^)]     (2.36) 

where  E°  is  the  standard  EMF  of  the  cell  and  k  is  again  writ- 
ten for  (RT  In  10)/F.   The  study  included  three  ionic  strengths, 
namely  0.3809,  0.6729,  and  0.8720  mol  kg-1  of  water,  corres- 
ponding to  seawater  with  salinities  of  20,  35,  and  45  /00 
respectively.   It  was  found  that  Earned 's  Rule  is  indeed 
valid  for  these  mixtures  in  the  temperature  range  5  to  50°C. 
By  casting  the  Nernst  equation  in  the  following  form: 
E  +  klogyHC1  =  E°  -  2klogIy°cl  +  2kayNaC1I       C2.37) 

=  A  +  3^NaCl 
the  Harned  slope,  a,  may  be  evaluated  by  the  method  of  least- 
squares.   It  is  then  possible  to  calculate  the  trace  activity 
coefficient  of  HC1  using  Earned' s  Rule: 

logyHC1  =  logy£cl  -  CayNaC1)I  (2.38) 

The  trace  activity  coefficient,  y_*  •    can  be  calculated  in 
this  manner  for  yHC1  when  yNaC1  =  1,  that  is,  when  NaCl  is  the 
only  solute;  it  should  approximate  the  yHC1  which  would  actu- 
ally be  measured  in  a  solution  of  pure  NaCl  at  ionic  strength, 

I,  containing  only  a  trace  of  EC1 .   These  calculated  values 

tr 
of  YHC1  at  tile  various  temperatures  and  ionic  strength  agree 

closely  with  the  "trace"  activity  coefficients  actually 

25 
measured  by  Khoo  et  al .    with  0.01m  EC1  in  artificial  seawater 
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without  sulfate,  although  they  are  somewhat  higher.   From  this 

it  may  be  deduced  that  the  decrease  in  Ytrn-i  from  Y£ri  is 

nearly  accounted  for,  in  seawater,  by  the  effects  of  NaCl 

alone. 

A  recipe  for  "sulfate  free"  seawater  suggested  by  Khoo 

is  as  follows: 

NaCl  =  0.46444m 
KC1  =  0.01058m 
CaCl9  =  0.01077m 
MgClg  =  0.05518m 

where  m  is  molality  and  the  ionic  strength  contribution  of 
Na2S04  has  been  replaced  by  NaCl.   Since  KC1  and  CaCl2  are 
rather  minor  components,  a  simpler  recipe  could  be  followed 
in  which  KC1  and  CaCl2  are  replaced  by  NaCl  and  MgCl2  respec- 
tively such  that: 

NaCl  =  0.47502m 
MgCl2  =  0.06595m 

The  total  ionic  strength  of  such  a  mixture  would  be  0.67287m 
(salinity  =  35  /00)  and  the  corresponding  ionic  strength 
fractions  are  yNaC1  =  0,706  and  7M  C1   =  0.294  such  that 

NaCl  •xMgCl9  =  1.  Now  for  a  three  component  system  such  as 
HCl-NaCl-MgCl2  in  water,  a  logical  extension  of  Earned' s  Rule 
might  be : 

logYHC1  =  logr|cl  -  aiyNaC1I  -  a2yMgC   I       (2.39) 

where  c^  is  the  Harned  slope  determined  by  the  Macaskill 
study  and  ct2  is  that  to  be  determined  from  a  similar  study  of 
the  HCl-MgCl2  system  in  a  cell  identical  to  2.35.   Equation 
2.39  assumes  a  linear  dependence  of  logy^^-.  on  the  ionic 
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strength  fractions  of  either  component.   Thus  (ytr  ) 

LTHCl;NaCl/MgCl2 

can  be  calculated  for  the  case  where  y„  _.  +  y,  „       =1  that 

^NaCl    yMgCl2    x'    za3-z 


tr 


logCYHC1)  =  logy°  -  [a1(0.706)  +  a2(0 .294)] I      (2.40) 

These  calculated  values  should  agree  closely  with  the  values 

t  r 
of  ^HCl^sw  measured  by  Khoo  if  our  assumptions  are  correct 

and  if  NaCl  and  MgCl2  are  the  major  influence  on  the  decrease 
°f  YHC1- 

It  might  also  be  fruitful  to  make  measurements  in  a  cell 
of  type  2.35  actually  containing  all  three  components  in 
water,  keeping  the  ratio  NaCl/MgCl2  constant  as  it  is  in  our 
recipe,  and  varying  only  (mNaC1  +  mMgC1  )  with  m^  at  con- 
stant ionic  strength.   In  this  case,  measurements  should  yield 

values  of  a  that,  when  put  in  equation  2.38  with  y  =  1,  give 

tr 
values  of  yhc1  identical  to  those  calculated  by  2.40.   This 

was  investigated  in  the  present  study.26 


CHAPTER  3 
CONCEPT  OF  pH 


Theoretical  Background 

In  the  preceding  chapter  we  have  seen  how  thermodyna- 
mics coupled  with  experimental  practice  can  provide  informa- 
tion concerning  mean  activities  and  activity  coefficients  of 
electrolytes.   However,  when  it  becomes  necessary  to  assign 

the  corresponding  values  for  individual  ions  one  faces  a 

27  28 
problem  for  which  there  is  no  direct  solution.   '    This 

problem  manifests  itself  very  clearly  in  th.e  determination 

9  29 
of  pa~  or  hydrogen  ion  activity.  ' 

Recalling  cell  2.7,  electromotive  force  measurements 

would  yield  E  in  terms  of  the  following: 

E  =  E°  -  k  log  aHmclTcl  C3.1) 

It  is  then  possible  to  write  an  expression  for  pa_.  =  -log  a„ 

in  the  following  manner: 

E    F° 
P^g  =  — 5 —  +  log  mcl  +  log  yci  C3.2) 

In  order  to  obtain  a  thermodynamically  sound  value  for  pa^, 

H 

a  knowledge  of  yc1  would  be  needed.   For  a  conventional  rep- 
resentation of  this  quantity,  the  Bates-Guggenheim  conven- 

.  .   30 

tion   was  proposed. 
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Because  of  the  success  of  the  Debye-Hiickel  equation  in 

fitting  activity  coefficient  data  at  low  ionic  strengths, 

Bates  and  Guggenheim  proposed  that  the  activity  coefficient 

of  chloride  ion  be  expressed  in  the  following  form  at  ionic 

strengths  up  to  0.1m: 

-A/T 

log  Yri  =  — C3,3) 

01    1+1.5/1 

The  choice  of  the  coefficient  1.5  in  the  denominator  is  rea^- 

sonable  in  that  it  assigns  y^,  nearly  the  same  value  as  the 

mean  activity  coefficient,  y+,    for  NaCl  up  to  0.1m,   This 

"convention"  has  been  instrumental  in  providing  a  basis  for 

the  NBS  standard  scale  for  the  pa-.  =  pHCS)  of  various  dilute 

buffer  substances.   If  equation  3,2  is  cast  as  follows: 

p^Vc^  =  T-^°  +  lQs  mci         C3-4) 

then  a  plot  of  this  acidity  function,  pCa„Yr-i),  as  a  function 

of  nip,  should  be  approximately  linear.   This  has  been  carried 

31  32 

out   '    for  a  variety  of  buffer  substances  to  which  chloride 

has  been  added,  and  at  the  point  of  intercept,  pCarjYp-,)0, 
equation  3.3  was  applied  to  obtain  values  of  pH(.S)  : 

pH(S)  =  paH  =  P(aHYcl)°  +  log  yc1  (-3.5) 

The  need  for  a  pH(S)  scale  is  apparent  from  the  stand- 
point that  it  is  often  desirable  to  evaluate  pH  changes  with 
precision  greater  than  0.1  unit.   An  internally  consistent 
scale  of  this  sort  enables  the  experimenter  to  calibrate  his 
pH  measuring  equipment  at  pH  values  both  above  and  below  that 
of  the  sample  by  proper  choice  of  NBS  standards. 
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Cells  With  Liquid-Junction 

In  practice,  most  pH  determinations  are  made  using  cells 
of  the  type: 

Reference  electrode  |  KCl(M)||Soln  X|Pt;  H2  or  glass  (.3.6) 

where  the  most  common  reference  electrode  is  the  calomel 
electrode  in  either  3 . 5M  or  saturated  KC1 .   The  double  ver- 
tical lines  indicate  a  liquid  junction  between  solution  X 
and  the  KC1  solution.   The  cell  is  first  calibrated  by  employ- 
ing a  pH  standard  buffer,  S,  in  place  of  solution  X.   For 

this  solution,  pH(S)  is  known  and  the  potential,  E  ,  between 

s 

the  hydrogen  sensing  and  the  reference  electrode  is  recorded. 

A  sample  solution  of  unknown  pH  is  introduced  into  the  cell 

replacing  the  standard  and  a  new  potential,  E  ,  is  measured. 

Ideally,  the  pH(X)  is  related  to  pH(S)  and  the  measured 

quantities  by  : 

(E   -  E  )F 
pH(X)  =  pHCS)  -   R*  ln  ^  (3.7) 

This  is  the  operational  definition  of  pH  and  provides  a 
simple  way  in  which  to  relate  the  pH  of  an  unknown  solution 
to  the  NBS  standard  scale.   Unfortunately,  pH(X)  values 
determined  in  this  manner  never  fall  exactly  on  the  conven- 
tional scale  because  of  a  residual  potential,  E.,  across  the 
liquid  junction  (vertical  lines  in  3.6).   This  potential  is 
a  complicated  function  of  the  activities  and  transference 
numbers  of  the  several  ionic  species,  i,  in  the  boundary 
layers  of  the  junction  and  may  be  represented  by: 
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2    t 
Ej  =  "f  {  ^dln  miYi  C3.8) 

where  t.  is  the  transference  number  for  species  i  and  z.  is 
its  charge.   An  exact  evaluation  of  E .  is  not  possible  with- 
out a  knowledge  of  a„  which,  in  turn,  is  the  object  of  the 
calculation.   This  is  a  dilemma  which  can  only  be  resolved 

from  "outside  the  realm  of  thermodynamics." 

33 
The  Henderson  method   of  integration  of  equation  3.8 

has  been  moderately  successful  in  the  estimation  of  the 
liquid-junction  potential.   Henderson  assumed  that  the 
junction  consists  of  a  continuous  series  of  mixtures  of 
solutions  1  and  2.   The  activity  of  each  ion  is  assumed  to 
be  equal  to  its  concentration,  and  its  mobility  is  constant 
throughout  the  junction.   The  concentration  of  each  species 
is  then  given  as  a  function  of  the  concentration  in  the  end 
solution  and  the  mixing  fraction,  a.   The  transference  num- 
ber is  given  in  terms  of  concentrations,  c,  and  mobilities, 
Ui,  as: 

c..u  . 

*i  =  ~ — —  C3.9) 

a  S  c . u .  +  (1  -  a)  Z   c.u. 
.   x  i  l  *■  l 

Substitution  of  3.9  into  3.8  and  integration  give: 

_  RT  z  cur^Hc':  -  c:>     z  clu. 

E-i  -  f n i In —        (3.10) 

If  the  ionic  mobilities  in  the  two  end  solutions  are  taken 
equal  to  the  mobilities  at  infinite  dilution  then  one  can 
write : 
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RT  cu1  -  v  )  -  (.u  -  v  )    u  +  v 

e.  =  S?  _i — 1 2 — 2_  ln  _i — i         C3>11) 

3  &i   ~   Vi>  "  <U2  "  V2}     U2  +  V2 

In  this  equation,  U  is  Z   c+^+  an(*  V  is  E  c_^°   f°r  the 
cations  and  anions  in  the  end  solutions  1  or  2 ;  U  and  V 
are  S  c+A°|z+|  and  Z   c_X°|z  |  respectively.   The  term  \°    is 
the  limiting  ionic  conductivity.   The  concentrations  are  in 
molarity.   Equation  3.11  often  yields  useful  estimates  of 
the  magnitude  and  sign  of  the  liquid-junction  potentials  and 
is  used  in  this  work. 

pH  in  Seawater 

Seawater  is  a  special  type  of  solvent  medium  in  that 
its  pH  is  regulated  within  narrow  limits.   Due  to  the  complete 
dissociation  of  the  primary  electrolyte,  NaCl ,  the  ionic 

strength  is  high  and  seawater  possesses  the  favorable  charae- 

34  35 
teristics  of  a  "constant  ionic  medium."   '    Media  of  this 

sort  can  be  expected  to  demonstrate  minimal  changes  in  Ytt 
as  small  changes  in  chemical  composition  occur.   Furthermore, 
if  pH  standards  were  available  which  match  the  ionic  strength 
of  seawater,  residual  liquid-junction  potentials  would  be 
expected  to  be  minimized.   Thus,  an  experimental  scale  of  pm^ 
could  be  set  up  rendering  the  possibility  of  measuring  hydro- 
gen ion  concentration  in  addition  to  activity.   Bates  and 
Mucaskill   have  shown  the  minor  variation  in  log  Yrrp-i  when 
NaCl  in  sulfate-free  seawater  is  replaced  by  small  quantities 
of  HC1  at  constant  ionic  strength.   In  addition,  it  appears 
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that  seawater  of  35  /00  salinity  indeed  effectively  nulli- 
fies the  residual  liquid- junction  potential. 

From  the  Nernst  equation  one  can  derive  the  following 
for  a  cell  of  type  2.7: 

E  =  E°  -  k  log  mHmcl  -  2k  log  yhc1         (.3.12) 

It  is  useful  to  alter  the  standard  state  such  that  the 

activity  coefficient  becomes  unity  at  zero  molality  of  H.C1 
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in  seawater  solvent  rather  than  in  pure  water.   '    This 

is  accomplished  through  the  following  substitutions  in 

equation  3.12: 

E°*  =  E°  -  2k  log  y^;  yEC1  -  Y^Y^        C3.13) 

tr 
where  YttC-,  is  the  transfer  activity  coefficient  or  medium 

effect  in  going  from  water  to  seawater  and  is  actually 

the  trace  activity  coefficient  of  HC1  in  seawater.   This 

yields : 

E  +  k  log  mHmcl  =  E°*  -  2k  log  YJC1        C3.14) 

The  quantity  E°*  may  be  determined  by  plotting  the  left  side 
of  equation  3.14  as  a  function  of  mFC1  when  small  quantities 
of  HC1  replace  NaCl  and  extrapolating  to  m  =  0.   These  types 
of  plots  have  been  found  to  be  straight  lines  with  the  slope 
amounting  to  only  0.04mV  at  0.01  mol  kg   of  added  SCI.   In 
view  of  this,  E°*  can  be  determined  routinely  with  0.01m  HC1 
with  a  reasonable  estimate  of  Y$C1  at  this  concentration. 
By  assuming  that  YgC1  is  virtually  constant  in  the  seawater 
medium  when  small  quantities  of  buffer  substances  are  added 
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in  place  of  NaCl ,  as  is  the  case  with  the  addition  of  HC1 , 
one  can  determine  pm„  =  -log  (nO  from  dilute  buffer  solu- 
tions in  seawater  by  measurements  of  cell  2.7  and  the 
relation: 

E  -  E°* 
pmH  =  — k —   +  log  mCl  (3.15) 

25 

Measurements  of  this  sort  have  been  carried  out   under  these 

conditions  and  point  to  the  practicality  of  a  standard  pm„(S) 
scale  in  seawater  from  which  values  of  prn^CX)  could  be  ob- 

n 

tained  through  use  of  the  operational  definition: 

(E   -  E  )F 
pm^X)  =  pmH(S)  -   R*  ln  fQ  (3.16) 

Since  nullification  of  the  residual  liquid-junction  potential 
requires  a  constant  ionic  medium  and  matched  ionic  strength 
between  the  standard  buffered  seawater  and  sample,  it  is  of 
interest  to  determine  the  magnitude  of  liquid-junction  errors 
due  to  salinity  variations.   This  is  investigated  in  the 
present  study  by  calibrating  a  cell  of  type  3.6  with  a  buff- 
ered seawater  at  35  /00  salinity  where  the  value  pm-.CS)  has 

n 
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been  determined.    Values  of  pmH(X)  can  be  calculated  by 

the  operational  definition  from  similar  measurements  above 
and  below  35  /00.  These  results  can  be  compared  with  pre- 
determined values  from  cells  without  liquid  junction. 


pH  in  Clinical  Media 

Measurement  of  pH  in  biological  fluids  is  an  important 

39  40 
part  of  clinical  diagnosis   '    and  is  normally  carried  out 
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with  a  glass  electrode  assembly  in  a  cell  such,  as  3.6. 
Standardization  of  the  cell  at  a  point  close  to  the  pH  of 
the  sample  is  possible  by  means  of  the  NBS  scale.   It  is 
also  desirable  to  match  the  ionic  strength  of  the  standard 
buffer  to  that  of  the  sample  to  minimize  residual  liquid- 
junction  effects.   As  was  previously  shown,  the  NBS  standard 
buffers  are  the  result  of  a  convention  proposed  for  use  at 
ionic  strengths  of  0.1m  or  below  and  for  that  reason  are  more 

dilute  than  is  desirable  for  clinical  work  where  the  ionic 

41 
strength  is  about  0.16m.  Attempts  have  been  made  to  es- 
tablish primary  reference  standards  at  this  ionic  strength 
using  a  solution  of  Tris-HCl  (0.05m),  Tris  (0. 01667m),  and 
NaCl  (0.11m);  however  pH(S)  values  assigned  to  this  buffer 
in  the  usual  manner  were  found  to  be  inconsistent  with  the 

NBS  scale  when  compared  in  cells  of  type  3.6.   Vega  and 
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Bates   have  proposed  the  use  of  two  substituted  ethanesul- 

fonic  acids  for  pH  control  in  the  physiological  range  of  pH 

and  at  I  =  0.16m  which  are  investigated  for  their  usefulness 
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in  the  present  work.  " 


CHAPTER  4 
OSMOTIC  COEFFICIENTS 


Theoretical  Background 

The  isopiestic  method  is  an  important  tool  for  studies 
of  electrolytes  in  solution  in  that  it  gives  a  direct  meas- 
ure of  the  osmotic  coefficient.   In  making  a  measurement,  a 
series  of  solutions  containing  the  salts  of  interest  are 
allowed  to  equilibrate  with  a  series  of  solutions  of  known 
vapor  pressure  at  the  vapor  pressure  of  water  through,  the 
vapor  phase.   The  activity  of  the  solvent  is  related  to  the 
vapor  pressure,  p,  by  the  following: 

aw  =  p/p°  C4.1) 

where  p°  is  the  partial  pressure  of  pure  solvent,  in  this 
case,  water.   Strictly  speaking,  this  relation  should  be  a 
ratio  of  fugacities;  however,  the  correction  is  nearly  iden- 
tical for  the  numerator  and  denominator  and  therefore  effec- 
tively cancels.   When  isopiestic  vapor  equilibrium  is  at- 
tained, a   is  the  same  for  all  solutions.   Recalling  equation 
1.9,  the  criterion  for  isopiestic  equilibrium  is; 

(vm*}test  =  (vm<j5)ref  (4'2) 

This  allows  one  to  calculate  d>     directly  if  A      „  is 

rtest        J  ref 

known.   For  this  reason,  a  reference  solution  for  which-  $ 
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is  known  is  always  equilibrated  with-  the  test  solutions. 
For  a  mixture  of  two  salts,  equation  4.2  is  simply: 

Cvlml  +  v2m2)c»test  =  (vm^ref  ^4.3) 

The  success  of  this  method  depends  on  three  essential  factors 

44 
outlined  by  Sinclair   and  mentioned  briefly  here: 

1)  Good  thermal  conduction  between  solutions  should 
be  provided  in  order  to  ensure  solvent  transfer 
rather  than  a  temperature  gradient. 

2)  Small  quantities  of  solution  nearly  at  isopiestic 
equilibrium  should  be  used  so  that  little  solvent 
need  actually  distill. 

3)  Stirring  should  be  sufficient  to  allow  for  maximum 
diffusion  of  solute  and  good  heat  conduction. 

As  mentioned  in  Chapter  1,  the  osmotic  coefficient  des- 
cribes the  non-ideal  behavior  of  the  solvent ,   This  non- 
ideality  is  largely  due  to  the  interaction  of  the  electrostat- 
ic  forces  of  the  ions.   Therefore,  as  a  first  approximation 
one  can  say  that  <j>  is  proportional  to  the  charge  density  of 
the  solution.   This  would  imply  that  two  solutions,  one  of 
NaCl  and  one  of  KC1 ,  at  the  same  ionic  strength,  would  have 
similar  osmotic  coefficients.   This  is  not  the  case,  however, 
as  a  1  molal  solution  of  NaCl  has  a  <J>  value  of  0.967  while 
at  the  same  concentration  of  KC1 ,  $  =  0.897.   In  this  sense, 
the  solvent  activity  must  be  controlled  by  other  more  specif- 
ic forces  dependent  on  the  nature  of  the  ions.   In  general, 
one  might  say  that  two  salts  of  the  same  charge  type  generate 
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different  osmotic  coefficients  as  a  consequence  of  specific- 
ities in  cation-anion  interactions  and/or  different  hydration 
characteristics . 

Solutions  of  Single  Electrolytes 

Since  the  work  of  Debye  and  Hxickel    there  have  been- 

45 
several  advances   in  explaining  theoretically  the  behavior 

of  electrolytes  at  appreciable  concentrations.   However, 

only  in  the  last  two  decades  has  extensive  correlation  of 

experimental  data  with  theory  resulted  in  useful  semi- 

.  -  18  46-^0 

empirical  treatments.   » 

From  equation  2.20  for  the  logarithm  of  the  activity  co- 
efficient according  to  Debye  and  Hiickel ,  we  have,  on  changing 

47  ^1 

the  symbolism  slightly,   ' 

S|z+z_|vT 

In  Y+  =  C4.4) 

1  +  A/I 

where  now  S  =  -1.17202  is  the  Debye  limiting  slope  at  25° C, 
/z+z_/  =  Z  is  the  charge  factor,  and  A  is  the  adjustable  term 
Ba.   If  this  relationship  is  differentiated,  substituted  into 
equation  1,19,  and  integrated  analytically,  the  following  ex- 
pression for  the  osmotic  coefficient,  $,    is  obtained: 

225*  -  l  +  8Z   [(1  +  A/I) 1 2  In  CI  +  A/I)]  C4.5) 

I        A°I  (1  +  A/T) 

For  higher  concentrations,  where  this  equation  would  normally 
fail,  additional  arbitrary  terms  may  be  added  in  higher  pow- 
ers of  ionic  strength,  such  that: 
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^L  =  !  +  _SZ  [C1  +  A/I)  _  1 2  In  (1  +  A/T)] 

I        AJI  (1  +  A/T) 

+  BI  +  CI2  +  DI3  +  EI4  (4.6) 

The  corresponding  equation  for  the  activity  coefficient  can 

be  written: 

ln-Y+  =  S +  (.2B)I  +  C3/2C)I^  +  C4/3D)I°  +  C5/4E)I^  C4.7) 

1  +  A/T 

Equation  4.6  can  be  applied  to  osmotic  coefficient  data  for 

solutions  of  single  electrolytes  and  the  parameters  A  through 

E  evaluated  by  the  method  of  least-squares.   The  parameters 

are  chosen  such  that  the  sum  squared  of  <f>  observed  minus  $ 

calculated  or  £  ($„ Ko  ~  ^o-i /•>)   is  minimized  over  the  ionic 

O  DS      Cct-L  C 

strength  range.   The  parameters  can  then  be  tabulated  for  a 
variety  of  electrolytes.   It  is  possible  to  employ  these  data 

in  Scatchard's  neutral-electrolyte  treatment  for  electrolyte 

50  52 
mixtures.   '    Although  this  method  represents  the  osmotic 

and  activity  coefficients  accurately  for  many  pure  electro- 
lyte systems,  there  is  little  promise  of  a  simple  physical 
interpretation  of  the  parameters  obtained  in  this  way. 

Pitzer,  on  the  other  hand,  has  taken  a  somewhat  differ- 
ent approach  in  seeking  to  use  fewer  parameters,  each  of 
which  retains  as  much  physical  meaning  as  possible.   By  mod- 
ification of  the  Debye-Hiickel  equation  in  the  manner  of 
Rasaiah  and  Friedman   '    and  using  a  three-term  power  series 
expansion,  he  derives  a  rigorous  expression  which  relates 
the  osmotic  coefficient  to  an  "electrostatic"  term  plus  a 


31 


series  of  virial  coefficients,  each,  of  which  is  a  peculiar 

45 
and  purely  theoretical  function  of  ionic  strength. 

In  these  terms,  the  properties  of  single  electrolyte 

solutions  take  the  following  forms  for  an  electrolyte  MX: 


-  1  =  |z+z_|f$  +  mC2vMvx/v)Bj 


MX 
+  m2[2(vMvx)3/2/v]C*x     C4.8) 

where         f*  =  — A  ^     _,  A*  =  -0.392 
1  +  1,2/1 

BMX  =  4x     +    SMX}  expC-2/1) 
It  can  be  seen  that  f   is  the  electrostatic  term  coming  di- 
rectly from  the  Debye-Hiickel  theory;  it  reduces  to  the  limit- 
ing slope  for  osmotic  coefficients  (-0.392)  at  low  concentra- 
tions.  The  parameter  B^L  represents  interactions  between 
ions  of  like  charge  and  may  be  represented  by  a  linear  term, 
SMX  '  plus  a  second  term,  8^  ^  exp(-2/T).   This  function  has 
the  property  of  a  rapid  change  with  /!  at  low  ionic  strengths 
and  a  smooth  approach  to  a  constant  value  at  high  ionic 
strengths.   The  third  parameter,  C*  ,  represents  triplet  in- 
teractions involving  two  ions  of  like  charge  and  one  of  oppo- 
site charge  and  is  certain  to  be  small  in  magnitude.   An  ex- 
pression analogous  to  equation  4.8  may  be  written  for  the 
activity  coefficient,  with  appropriate  changes  in  these  virial 
coefficients : 

In  Y±  =  |z+z_|fY  +  m(2vMvx/v)B^x  +  m2 [ 2(vMvx)3/2/v] C^  C4.9) 
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where 
fY   =      A* [/I /CI   +   1.2/1)   +    C2/1.2)    In    (1   +   1.2/T)] 

BMX  -   BMX  +   VI   /   B^Cx)dx 

^Mx        J/^MX 

Thus,  osmotic  coefficient  data  fitted  to  equation  4.8  may  be 
converted  readily  into  activity  coefficient  data  for  the  sol- 
ute. 

Solutions  of  Mixed  Electrolytes 

In  seeking  to  describe  the  variation  of  osmotic  coeffi- 
cients of  mixtures  of  electrolytes,  two  noteworthy  treat- 
ments have  been  used  extensively,  namely  Scatchard's  neutral- 

50  52  55—57 

electrolyte   '   '      treatment  and  the  equations  of 

Pitzer.46'58"60 

As  mentioned  above,  Scatchard's  equations  for  mixtures 
of  electrolytes  result  directly  from  single  electrolyte  equa- 
tions 4.6  and  4.7  where  the  Debye-Hiickel  theory  was  extended 
using  power-series  terms  in  ionic  strength.   Rearrangement  of 
Scatchard's  general  equation  for  the  osmotic  coefficient  of  a 
mixture  of  two  1:1  electrolytes  in  water  yields  the  following: 

2*  =  2*Pa  +  2*&b  +  yAyBia>oi  +  b02x  +  bo3l2> 

+  CyA  -  yB)yAyBKb12i  +  b13i2)  C4.10) 

where   <f>°  is  the  osmotic  coefficient  in  the  pure  single  elec- 
trolyte solution  derived  from  equation  4.6,  y.  is  the  ionic 
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strength,  fraction  of  component   j,  and  th.e  terms  b..  are 
various  mixing  parameters.   Osmotic  coefficients  obtained 
from  such  mixtures  in  conjunction  with  the  best -fit  param- 
eters from  equation  4.6  may  be  used  to  evaluate  the  mixing 
parameters,  b. .,  by  the  method  of  least -squares.   Correspon- 
ding equations  can  also  be  written  for  the  activity  coeffi- 
cients of  the  individual  species,  which  may  be  calculated 
from  the  appropriate  forms  when  these  mixing  parameters  are 
available.   Similarly,  one  can  express  the  osmotic  coefficient 
for  a  mixture  of  1:1  and  1:2  electrolytes  in  water  by: 

(1  +  yA)<t>  =  2cf)°y4  +  (ji  y   +  mixing  terms       C4.ll) 

Pitzer's  general  equation  for  electrolyte  mixtures  is 
not  readily  reduced  to  terms  in  <j>°  and  y  for  the  single 
electrolytes.   For  a  solution  composed  of  two  electrolytes 
MX  and  NX  with,  a  common  ion  X,  with  no  restriction  on  charge, 
his  equation  rather  takes  the  form: 


4  -  1  =  [2/Cffljj  +  nijj  +  mx)3  {If*  +  mMmxB|x  +  ny^R 
+  VN^MN  +  I0MN)  +  •  VZM|1/2mMmXC!? 


NX 


+   z 


MX 
l/2m   2, 


x/zN|    m^xCjx  +  nyn^^}  C4.12) 

where  0j_  and  i „NX  arise  from  a  consideration  of  additional 

interactions  among  combinations  of  ions.   The  term  0,' 

MN 

arises  from  the  dependence  of  0„N  on  ionic  strength  and  is 
expected  to  be  small.   The  other  terms  are  identical  to  those 
in  equation  4.8  and  are  determined  from  data  on  the  single 
electrolyte  solutions,  MX  and  NX. 


CHAPTER  5 
EXPERIMENTAL 


Electromotive  Force  Me  astir  ements 

Materials 

Purified  hydrochloric  acid  for  the  EME  measurements  was 
obtained  by  distilling  reagent  grade  CFisher)  HC1 ,  retaining 
the  middle  fraction  from  each  of  two  distillations.   The 

product  was  analyzed  gravimetrically  for  chloride  by  precip- 
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it at ion  of  AgCl  using  a  standard  procedure.     Results  of 

analysis  indicated  5.63600  mol  chloride  per  kg  Cair  weight) 
of  stock  solution.   A  secondary  stock  (0.186625  mol  kg-1  sol- 
ution) was  prepared  by  dilution  and  this  was  used  in  prepara- 
tion of  the  mixtures. 

A  solution  of  magnesium  chloride  was  prepared  by  dissolv- 
ing  reagent  grade  (Mallinckrodt )  MgCl2  •  6H20  in  water  to 
make  a  solution  approximately  1.5  molal.   The  resulting  solu- 
tion was  analyzed  gravimetrically  to  yield  a  stock  solution 
3.04958  mol  kg"1  in  chloride. 

Sodium  chloride  was  purified  by  twice  recrystallizing 
reagent  grade  (Mallinckdrodt )  NaCl  from  water  and  drying  the 
product  at  300°C  prior  to  use.   Stock  solutions  were  prepared 
by  weighing  the  dried  solid. 
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Phosphate  buffer  solutions  for  the  pH  measurements  were 
prepared  from  the  NBS  Standard  Reference  Materials:   KHgPO 
SRM186Ic,  and  Na2HP04 ,  SRM186IIc,   The  solids  were  dried 
overnight  at  110°C  prior  to  use  and  the  solutions  prepared 

by  weighing  the  dried  solids  and  adding  degassed  water 

9 
according  to  standard  procedure. 

Buffers  of  N-tris(hydroxymethyl)methyl-2-aminoethane- 
sulfonic  acid  (TES)  and  a  N-2-hydroxyethylpiperazine-N' -2- 
ethanesulfonic  acid  (HEPES)  were  prepared  from  the  commercial 
solids  (Sigma)  which  had  been  recrystallized  twice  from  80% 
ethanol-water  and  dried  at  80°C  in  a  vacuum  desiccator,   Each. 
was  then  weighed  accurately  into  a  glassy-stoppered  flask  to 
which  was  added  an  appropriate  amount  of  standardized  NaOH 
solution.   Dried  NaCl  was  then  added  as  necessary  to  achieve 
the  desired  ionic  strength. 

Standard  carbonate-free  sodium  hydroxide  solution  was 
prepared  by  dilution  of  a  50%  (wt/wt)  NaOH  solution  with  de- 
gassed water.   It  was  standardized  by  titrating  with  primary 
standard  grade  (NBS  SRM84h)  potassium  acid  phthalate. 

A  buffer  solution  of  tris(hydroxymethyl)aminomethane 
(Tris)  was  prepared  from  a  solution  of  the  primary  standard 
grade  (Sigma)  base  to  which  was  added  standard  HC1  solution 
and  the  appropriate  amount  of  NaCl . 

Potassium  chloride  for  the  synthetic  seawater  solutions 
was  reagent  grade  (Mallinckrodt )  which  had  been  recrystal- 
lized twice  from  water  and  dried  overnight  at  120°C,   Solu- 
tions were  prepared  by  weighing  the  dried  solid. 
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A  stock  solution  of  calcium  chloride  was  prepared  from 
reagent  grade  CMatheson ,  Coleman,  and  Bell)  CaCl2  •  2H„0  by 
dissolving  the  salt  in  water  to  make  a  solution  approximately 
2.0m.   The  solution  was  analyzed  gravimetrically ,  and  the 
results  indicated  2.73604  mol  kg    solution. 

Sodium  sulfate  was  reagent  grade  (Fisher)  anhydrous 
Na2S04  dried  at  120°C  in  vacuo  for  at  least  6  hours.   A 
stock  solution  was  prepared  by  weighing  the  dried  solid. 

Water  used  in  the  preparation  of  all  solutions  was  de- 
ionized  in  a  central  building  source  and  redistilled  in  our 
laboratory  from  an  all-glass  still.   For  buffer  solutions 
the  water  was  degassed  either  by  boiling  or  bubbling  N„ . 

Hydrogen  gas  was  obtained  in  a  commercial  cylinder  and 
purified  by  passage  over  a  platinum  "DeOxo"  catalyst  prior 
to  introduction  into  the  cells. 

Silver  oxide  for  the  preparation  of  the  silver-silver 
chloride  electrodes  was  prepared  by  addition  of  a  solution  of 

sodium  hydroxide  to  silver  nitrate  solution  as  described  by 

9 
Bates.    The  precipitate  was  washed  30  times  with  distilled 

water , 

Platinizing  solution  for  the  hydrogen  electrodes  was 
prepared  by  dissolving  2  grams  of  hexachloroplatinic  acid  in 
100ml  of  2M  HC1. 

Mercury  for  the  calomel  electrode  was  commercial  in- 
strument mercury  (Bethlehem)  triply  distilled  in  continuous 
vacuum.   The  calomel  was  reagent  grade  (Fisher)  Hg0Cl 
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Preparation  of  Cell  Solutions 

All  solution  mixtures  were  prepared  by  weighing  and 
combining  stock  solutions  containing  a  known  amount  of  each- 
substance  such  that  the  desired  mole  ratios  were  achieved. 
The  stock  solutions  were  made  sufficiently  concentrated  so 
that  water  was  always  added  as  a  last  step  to  reach  the 
desired  ionic  strength.   Weighings  were  performed  on  a 
Sartorius  5-decimal  or  a  Mettler  3-  or  4-decimal  balance  as 
appropriate  such  that  accuracy  was  always  better  than  1  part 
in  ten  thousand.   Vacuum  corrections  were  made  in  all  weigh- 
ings. 

For  the  HCl-MgCl„  mixtures,  HC1  was  weighed  into  a  tared 
glass-stoppered  erlenmeyer  flask  and  the  amount  of  MgCl„ 
needed  to  achieve  approximate  mole  fractions  of  0.9,  0.7, 
0.5,  0.3,  and  0.1  in  MgCl„  was  added.   One  solution  of  pure 
HC1  was  prepared  at  each  of  the  ionic  strengths.   The  quanti- 
ties of  HC1  stock  weighed  out  initially  were  such  that  the 
final  solutions  would  weigh  about  100  grams.   This  procedure 
was  carried  out  to  give  solutions  at  four  ionic  strengths, 
namely  0.1000,  0.3809,  0.6729,  and  0.8720  mol  kg"1. 

Three-component  mixtures  of  HCl-NaCl-MgCl2  were  prepared 
similarly  from  HC1  stock  and  a  solution  of  NaCl  and  MgCl„  of 
a  fixed  molal  ratio,  NaCl/MgCl2  =  7.202.   Four  solutions  were 
prepared  in  which  NaCl  and  MgCl2  together  contributed  frac- 
tions of  approximately  0.5,  0.3,  0.1,  and  0.0  to  the  total 
ionic  strength,  I  =  0.6729  mol  kg"  . 
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The  compositions  of  the  buffer  solutions  used  for  pH 
measurements  are  given  in  Table  5.   In  cases  where  buffer 
and/or  the  NaCl  concentration  was  varied,  the  buffer  was 
appropriately  diluted  and  NaCl  added  by  weighing  the  dried 
solid. 

Synthetic  seawater  solutions  buffered  with,  an  equimolal 
Tris  buffer  were  prepared  at  salinities  of  4Q,  35,  and  3Q°/QQ 
with,  three  buffer  concentrations,  0.06,  0,04,  and  0,02m,  at 
each,  salinity.   This  was  achieved  by  combination  of  a  concen- 
trated seawater  stock  solution  containing  less  than  the  de- 
sired amount  of  NaCl,  a  stock  solution  of  NaCl,  and  an  equi- 
molal Tris:Tris-HCl  buffer.   The  synthetic  seawater  stock 
solution  was  prepared  according  to  the  recipe  given  by  Khoo 
et_al. ,  such  that  if  all  of  the  required  NaCl  had  been  added 
the  salinity  would  have  been  greater  than  45°/00.   The  solu- 
tions were  prepared  by  weighing  out  an  aliquot  of  this  sea- 
water stock  and  calculating  the  amount  of  water  necessary  to 
dilute  it  to  the  desired  salinity.   Appropriate  quantities 
of  Tris:Tris-HCl  and  NaCl  solutions  were  then  added  by  weight 
and  finally  the  mixture  diluted  with  water.   It  should  be 
noted  that  Tris :Tris-HCl  was  added  at  the  expense  of  NaCl  to 
maintain  the  proper  chloride  content. 

Preparation  of  Electrodes 

The  hydrogen  electrodes  were  made  from  platinum  foil 
approximately  1cm2  in  area,  spot -welded  to  a  platinum  wire 
sealed  in  glass  tubing,   The  foil  was  platinized  in  a 
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solution  of  2%  Cwt/vol)  chloroplat inic  acid  in  2M  HC1  by 
electrolyzing  at  about  30mA  for  5  minutes.   Care  was  taken 
such  that  no  bubbles  formed  at  the  electrode  surface  during 
electrolysis.   The  electrodes  were  stored  in  distilled  water 
and  cleaned  after  use  by  boiling  in  "aqua  regia." 

The  silver-silver  chloride  electrodes  were  of  the 

9 
thermal-electrolytic  type   and  were  prepared  by  thermal 

decomposition  of  a  paste  of  silver  oxide  and  water  formed 
at  the  base  of  the  electrode.   The  base  was  a  helix  of  No,  26 
platinum  wire  about  2mm  in  diameter  sealed  in  glass  tubing. 
The  surface  of  the  silver  formed  was  then  converted  to  silver 
chloride  by  electrolysis  in  1M  HC1 .   The  electrodes-  were 
stored  in  dilute  (0.1m)  NaCI  and,  after  use,  were  cleaned  by 
boiling  in  concentrated  nitric  acid.   Following  preparation, 
the  electrodes  were  checked  for  bias  by  measuring  the  EMF  of 
each  electrode  versus  an  arbitrary  reference  electrode,   Any 
bias  was  taken  into  account  in  the  measurements. 

The  calomel  electrode  in  the  cell  with  liquid  junction 
was  prepared  with  3.5M  KC1 .   A  mercury  pool  about  1cm  deep 
was  formed  in  the  right-hand  half  cell  and  was  covered  with 
a  2mm  layer  of  calomel  which  had  been  wetted  with,  the  KC1 
solution.   The  cell  compartment  was  carefully  filled  to  the 
14/20  ground  glass  joint  at  the  top  and  a  platinum  wire 
sealed  in  glass  tubing  with  a  14/20  adaptor  was  carefully 
lowered  through  the  KC1  solution  and  calomel  into  the  mer- 
cury to  avoid  mixing  of  the  layers.   The  cell  compartment 
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was  completely  filled  with  KC1  solution  prior  to  this  step 
to  exclude  all  air  and  create  an  airtight  seal  at  the  joint. 

Equipment  and  Procedures 

The  cell  without  liquid  junction  is  shown  in  Figure  1 
and  is  described  by  Gary,  Bates,  and  Robinson.62   It  is  an 
all-pyrex  cell  except  for  teflon  stopcocks.   A  three-stage 
saturation  process  takes  place  in  the  compartments  to  the 
right  of  the  electrodes  and  insures  that  no  change  in  solu- 
tion composition  will  occur  in  the  electrode  compartments 
due  to  the  passage  of  hydrogen  gas. 

The  cells  were  first  flushed  with  hydrogen  to  remove 
any  traces  of  oxygen.   The  solutions,  also  deaerated  and 
saturated  with  hydrogen,  were  then  introduced  into  the  cell 
via  the  y-tube,  care  taken  so  as   not  to  introduce  oxygen. 
The  electrodes  were  then  placed  into  their  respective  com- 
partments.  The  cells  were  supported  in  brass  frames  and 
placed  in  a  water  bath;  temperature  control  to  ±0.01°C  in 
the  range  5  -  50°C  was  provided.   Purified  R"2  was  allowed 
to  bubble  constantly  through  the  saturator  over  the  hydro- 
gen electrode  and  out  of  the  gas  trap.   The  measurements 
were  shown  to  be  consistent  after  2  hours  at  25°C,  at  which 
time  the  temperature  was  lowered.   Equilibrium  was  reattained 
rapidly,  within  30  minutes,  and  the  next  measurement  re- 
corded.  This  procedure  was  followed  at  5°  intervals  to  5°C, 
and  then  the  temperature  was  cycled  upward  to  25°C  where 


Figure  1.   EMF  cell  without  liquid  junction 


42 


43 

the  original  measurement  was  checked  for  reproducibility. 
Measurements  were  then  made  at  5°  intervals  to  50°C  and 
back  to  25°C  for  a  final  reading.   Cell  measurements  ob- 
tained in  this  manner  demonstrated  excellent  reproducibility, 
within  0.04mV,  when  EMF's  were  corrected  to  standard 
1  atmosphere  pressure. 

The  cell  with  liquid  junction  is  shown  in  Figure  2.   A 
calomel  electrode  was  prepared  in  the  right-hand  compartment, 
and  the  KC1  solution  was  constantly  renewed.   The  left-hand 
compartment  was  equipped  with,  a  gas  tube  through  which  pre- 
saturated  H„  was  bubbled  over  the  electrode.   The  electrode 
was  seated  in  a  14/20  ground-glass  joint.   Gas  was  allowed 
to  escape  through  the  gas  trap.   Solutions  that  had  been 
deaerated  and  saturated  with  H»  were  introduced  into  the 
left-hand  compartment  followed  by  the  electrode,  care  being 
taken  to  avoid  air.   The  KC1  solution  being  the  heavier, 
formation  of  the  liquid-junction  above  the  3-way  stopcock 
was  permitted.   A  water  jacket  completely  surrounded  the 
cell  compartments,  and  each  solution  was  allowed  to  ther- 
mally equilibrate  for  1  to  2  hours.   After  formation  of  the 
junction,  the  calomel  and  hydrogen  electrodes  were  connected 
to  a  digital  voltmeter  and  the  EMF  recorded  when  it  appeared 
stable,  usually  within  5  minutes.   Inasmuch  as  the  EMF  was 
not  always  reproducible  from  day  to  day,  all  measurements 
were  normalized  to  a  fixed  value  for  a  phosphate  reference 
buffer,  measurements  of  which  were  constantly  repeated  through- 
out the  course  of  the  experiments. 


Figure  2.   EMF  cell  with  liquid  junction 
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Isopiestic  Measurements 

Materials 

Strontium  chloride  was  purified  by  twice  recrystalliz- 
ing  reagent  grade  (Mallinckdrodt \   SrCl2  from  slightly 
acidified  water.  A  stock  solution  was  prepared  and  analyzed 
gravimetrically  for  chloride.   The  results  indicated  1.63068 
mol  SrCl„  per  kg  solution. 

Sodium  carbonate  (.Fisher)  was  heated  at  300 °C  for  24 
hours  to  drive  off  H„0  and  was  then  dissolved  in  degassed 
water  to  yield  1.65165  mol  kg"   solution. 

Sodium  chloride  for  use  as  the  isopiestic  reference 
was  recrystallized  twice  from  water  and  dried  at  300°C. 
Solutions  were  made  by  weighing  the  dried  solid. 

Preparation  of  Solutions 

All  solution  mixtures  were  prepared  by  combining  stock 
solutions  by  weight  to  achieve  the  desired  mole  ratios. 
The  water  was  previously  degassed  to  minimize  sputtering  of 
the  solutions  when  they  were  subjected  to  vacuum.   The 
NaCl-SrCl2  and  the  NaCl-Na2C03  mixtures  were  prepared  by 
identical  methods,  with  the  ionic  strength  fraction  of  the 
second  component,  y„,  varying  from  0.17  to  0.85.   The  solu- 
tion with  the  highest  yR  also  necessarily  had  the  highest 
ionic  strength.   The  exact  compositions  of  the  solutions 
are  given  in  Tables  10  and  13. 
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Equipment  and  Procedures 

The  apparatus  for  the  isopiestic  measurements  consists 
of  a  vacuum  desiccator  containing  up  to  12  gold-plated  sil- 
ver cups.   The  cups,  containing  2  to  3  grams  of  solution, 
were  seated  on  a  copper  block  2.5cm  thick  to  ensure  thermal 
equilibration.   The  desiccator  was  slowly  evacuated  to  a 
pressure  of  about  25  torr,  carefully  avoiding  sputtering  of 
the  solutions.   The  entire  desiccator  was  then  immersed  in 
a  water  bath  maintained  at  25  ±  0.01°C  and  was  rocked  gently 
by  mechanical  means.   Stirring  was  accomplished  by  the  use 
of  either  glass  beads  or  platinum  gauze  in  each  cup.   Con- 
centrated solutions  (>lm).  required  only  3  days  to  a  week  to 
reach  equilibrium,  while  more  dilute  solutions  would  require 
up  to  2  weeks.   Duplicate  cups  were  run,  that  is,  2  cups  for 
each  solution.   At  the  end  of  the  equilibrating  period,  the 
cups  were  weighed  on  a  microbalance.   Knowing  the  weight  of 
the  solid  in  each  cup  and  the  weight  of  the  cup,  the  initial 
calculation  was  a  measure  of  wt .  solid/wt.  solution  from 
which  the  molality  of  each  solution  was  determined.   Equilib- 
rium was  assumed  complete  when  this  result  agreed  within 
0.1%  between  2  cups  containing  identical  solutions.   Follow- 
ing the  weighings,  the  cups  were  returned  to  the  desiccator 
and  the  amount  of  water  necessary  to  dilute  the  solutions 
to  a  lower  ionic  strength  was  calculated  and  added  dropwise 
with  a  graduated  syringe.   The  system  was  then  re-evacuated 
and  the  procedure  repeated  until  the  desired  ionic  strength 
range  was  covered . 


CHAPTER  6 
RESULTS 


Electromotive  Force  Measurements 


Standard  Potential  of  the  Silver-Silver  Chloride  Electrode 

Before  analyzing  results  from  EMF  Measurements  in  a 
cell  of  type  2.7,  one  needs  to  evaluate  the  quantity  E°  in 
the  Nernst  equation,  that  is,  the  standard  potential  of  the 
silver-silver  chloride  electrode.   It  has  been  discovered 
that  values  of  EV.AffC1  determined  in  different  laboratories 
lead  to  varying  results,  presumably  due  to  variations  in 
preparative  techniques.    For  this  reason,  routine  standard- 
ization should  be  continued,  and  such  determinations  have 
been  made  in  the  present  work. 

From  the  Nernst  equation  describing  the  EMF  of  cell 
2.7,  one  can  write: 

-  _  ™    4.605RT  ,.,  .        .        ,,.-,, 

E  -  E° jp—  (log  mHC1  -  log  yhc1)       (6.1) 

As  pointed  out  in  Chapter  2,  log  yhc1  can  be  accurately  rep- 
resented at  low  concentrations  by  an  equation  of  the  form: 

Am1/2 
"log  YHC1  =  TV"  1/2  "  3m  C6.2) 

1  +  m  ' 

where  A  is  the  constant  of  the  Debye-Hiickel  limiting  law  for 

a  particular  temperature.   Thus,  from  equation  6.1,  one  has: 
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AmV2 
1  +  m1^ 


E  +  2k  log*--.  --^-1/Q  U  E°  -  Sm 


HC1  ~  I 1/2  /=  E°  -  gm        C6.3) 


A  plot  of  the  left  side  of  equation  6.3  as  a  function  of  m 
should  give  a  straight  line  with  an  intercept  of  E° .   Bates 
and  Bower  "  have  tabulated  values  of  E°  at  various  tempera- 
tures based  on  extrapolations  of  this  sort,  only  with  the  in- 
clusion of  the  ion-size  parameter  which  can  be  varied  to  give 
the  best  fit  of  the  data.   Having  obtained  the  best  values 
for  E°,  one  can  then  return  to  equation  6,1  and  calculate 
values  of  log  yHC1 .   Considering  cell  2,7  containing  0.01 
mol  kg   HC1,  one  may  express  the  standard  potential  as: 

E°  -  E  +  2k  log  (0.01)yHC1  (6-4) 

Thus,  E°  may  be  calculated  over  a  range  of  temperatures  from 
one  measurement  only  at  each  temperature,  with  a  knowledge 
of  Yrci  at  that  concentration.   Table  1  gives  the  values  of 
E  and  E°  obtained  from  measurements  of  this  sort  in  the 
temperature  range  5  to  45°C.   Values  of  E°  from  Bates  and 
Bower's  extrapolations  are  given  for  comparison. 

Activity  Coefficient  of  HC1  in  Seawater 

In  an  attempt  to  elucidate  activity  coefficient  behavior 
of  HC1  in  saline  mixtures  approximating  seawater,65'66  EMF 
measurements  have  been  carried  out  in  mixtures  of  HCl-MgCl 
and  HCl-NaCl-MgCl2.   The  sodium  and  magnesium  chlorides,  the 
major  salt  components  of  seawater,67  exert  a  substantial 
influence  on  the  observed  behavior. 


Table  1.   Values  of  E  and  E°  for  th,e  Cell: 
Pt;  H2|HCl(0.01m)|AgCl;  Ag 
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t(.°C) 

E  (V) 

E°  CV) 

E°  (V) 

5 

+0,45951 

+0.23410 

+0,23413 

10 

0.46091 

0.23141 

0.23142 

15 

0.46216 

0.22857 

0,22857 

20 

0.46325 

0.22556 

0.22557 

25 

0.46417 

0.22238 

0.22234 

30 

0.46498 

0.21907 

0 . 21904 

35 

0,46565 

0,21566 

0.21565 

40 

0.46615 

0.21208 

0.21208 

45 

0.46665 

0.20836 

Q. 20835 

64 
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The  measurements  were  made  in  cells  of  the  type: 

Pt;  H2(.g,  latm>|HClCmA)  +  MgCl2CmB)  |  AgCl ;  Ag    (.6.5) 

and 

Pt;  H2Cg,  latm)|HClCmA)_  +  [NaCl/MgClgJ  Cnig)  |  AgCl ;  Ag    (.6.6) 

where  mA  and  mg  are  respectively  the  molality  of  HC1  and  the 
molality  of  added  MgCl2  or  the  mixture  NaCl/MgCU.   In  the 
first  case,  the  ratio  of  acid  to  salt  was  varied  while  the 
total  ionic  strength  was  kept  constant.   The  measurements 
were  carried  out  in  this  manner  at  four  ionic  strengths, 
namely  0.1,  0.3809,  0.6729,  and  0.8720  mol  kg-1.   The  three 
higher  ionic  strengths  correspond  to  seawater  of  salinities 
20,  35,  and  40  /00.   In  the  second  series  of  measurements, 
the  ternary  mixture  HCl-NaCl-MgCl2  was  investigated,  where 
the  ratio  of  HC1  to  total  salt  was  varied  at  an  ionic 
strength  of  I  =  0.6729  mol  kg-1  CS  =  35°/00) .   The  molal 
ratio  of  NaCl/MgCl2  was  always  kept  constant  at  7.202,  which 
would  be  the  same  ratio  as  in  seawater  consisting  of  only 
those  two  salts,   All  measurements  were  made  at  nine  temper- 
atures, 5  to  45  C.   The  results  of  all  of  the  EMP  measure- 
ments are  given  in  Table  2. 

The  HCl-MgCl2  mixtures  have  a  total  ionic  strength 
given  by  I  =  mA  +  3mB,   The  ionic  strength  fraction,  y   of 
MgCl2  is: 

yB  =  3mB/I  C6.7) 
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Thus,  both  m.  and  m„  can  be  calculated  from  the  values  of 
I  and  yB  given  in  Table  2.   For  the  ternary  mixtures  where 
mNaCl  =  7*202mM2:Cl  '  the  to"ta-'-  i°nic  strength  is  given  by: 

I   =  mA  +  1.2348mB  (6.8) 

where  yA  =  1.2348mB/I  (6.9) 

Mean  ionic  activity  coefficients  of  HC1 ,  y.,    in  the  pres- 
ence of  MgCl2  are  easily  calculated  from  the  Nernst  equation 
expressed  in  terms  of  y_. : 

log  YA  =  E°2k  E  "  °'5  log  (1  "  yB)(1  "  yB/3)l2     (6-10) 
and  for  the  ternary  mixtures: 

l0g  YA  =  E°2k  E  "  °-5  log  (1  _  yB)(1  ~  °-09802yB)I2    (6.11) 

Just  as  in  equation  2.37,  the  Harned  slope  may  be  calculated 

by  the  method  of  least-squares  for  any  given  values  of  I  and 

yB.   With  modification  for  HCl-MgCl2,  the  equation  takes  the 

form: 

E   +   k   log    (1    -   yB)(l    -   yB/3)    =   E°    -    2k    log    Iy°4   +    2kctyBI 

=    A   +   ByB  (6.12) 

where  A  and  B  are  determined  by  linear  regression  analysis. 
For  the  ternary  mixture,  equation  6.12  takes  the  form: 

E  +  k  log  CI  -  yB)d  -  0.09802yB)  =  A  +  ByR      (.6.13) 

The  results  for  a  (=B/2kI)  and  log  y*    are  given  in  Table  3, 
along  with  the  standard  deviations  of  fit  at  each  temperature 
and  ionic  strength.   Values  of  the  temperature  variation  of 
a,  da/dT,  are  given  at  the  foot  of  Table  3. 
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pH  Standardization  in.  Seawater 

In  an  attempt  to  estimate  the  magnitude  of  liquid- 
junction  errors  inherent  in  pH  measurement  in  seawater  of 
varying  salinity,  measurements  in  a  cell  with  liquid  junction 
have  been  carried  out  in  buffered  artificial  seawater  at 
salinities  of  30,  35,  and  45°/0O,  each  containing  Tris  buffer 
at  molalities  of  0,02,  0,04,  and  0,Q6. 

The  results  of  these  measurements  were  treated  in  two 
ways,  using  the  operational  definition  of  pH  Cequation  3.161. 
First,  values  of  pn^CXl  were  calculated  where  E  and  pm_CSl 
referred  to  the  35°/00  =  S,  0.04m  Tris  solution.   Values  of 
pmHCS)  were  obtained  from  Ramette  et  al . 34   In  the  second 
case,  pH^gg  was  calculated  by  referencing  E  to  the  E   and 
pH(.S)  values  for  the  NBS  1:3.5  phosphate  buffer.   The  results 
from  these  measurements  and  calculations  are  given  in  Table  4. 

pH  Standardization  in  Clinical  Media 

Attempts  to  establish  a  primary  reference  standard  for 
the  measurement  of  pH  near  7.4  at  an  ionic  strength  of  0.16, 
corresponding  to  blood  plasma  and  other  clinical  media, 
have  met  with  difficulty.   Ideally,  one  would  like  to  have 
a  standard  with  pH  at  or  near  that  of  the  sample.   Further- 
more, the  ionic  strength  of  the  standard  should  match  that 
of  the  sample.   Although  buffers  have  been  proposed  which 
meet  these  requirements,41'68  they  have  been  shown  to  be  in- 
consistent with  the  NBS  primary  scale,  that  is,  the  pH 
values  assigned  based  on  measurements  in  cells  without  liquid 
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Table  4.   Comparison  of  pmH(X)  and  pHNBg  Obtained  by  EMF 
Measurements  of  the  Cell  Hg;  Hg2Cl2 ,  3.5m  KCl||Soln. 
S  or  X  |  H2  (g,  latm);  Pt  at  25° C  With  Corresponding 
pmH(S)  Values. 

Soln.  X  or  S  S(°/00)  -E(V)  pnVgOO  pH^g  pn^CS)34 

sw  +  0.02m  Tris    30  0.73473  8.190  8.198  8.185 

sw  +  0.04m  Tris    30  0.73506  8.195  8.203  8.185 

sw  +  0.06m  Tris     30  0.73494  8.193  8.202  8.187 

sw  +  0.02m  Tris     35  0.73515  8.197  8.205  8.198 

sw  +  0.04m  Tris    35  0.73538  8.201  8.209  8.201 

sw  +  0.06m  Tris    35  0.73522  8.198  8.206  8.201 

sw  +  0.02m  Tris     40  0.73552  8.203  8.211  8.216 

sw  +  0.04m  Tris    40  0.73563  8.205  8.213  8.217 

sw  +  0.06m  Tris     40  0.73563  8.205  8.213  8.220 

1:3.5  Phosphate  I  =  0.1m  0.68841  (pH(S)  =  7.415) 
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junction  do  not  agree  with,  the  pH.  values  determined  by- 
measurements  in  cells  with  liquid  junction  and  the  opera- 
tional definition  (equation  3.7).   Pursuant  to  this  problem, 
extensive  measurements  involving  the  buffers  of  TES,  HEPES , 
and  Tris  have  been  carried  out  in  conjunction  with  the  NBS 

1:3,5  phosphate  "blood  buffer"  in  an  effort  to  reveal  the 

43 
causes  of  the  discrepancies.     Factors  investigated  were 

a)  the  possible  failure  of  the  convention  for  log  Yp.  at 
ionic  strengths  above  0,1m,  and  b).  the  appreciable  residual 
liquid-junction  potential  when  the  relatively  dilute  NBS 
reference  solutions  are  replaced  by  buffers  at  I  =  0.16m. 
Measurements  have  been  made  in  cells  without  liquid 
junction  to  determine  the  EMF  for  solutions  of  the  NBS 
phosphate  "blood  buffer"  containing  0.01,  0.02,  0.04,  and 
0.06  mol  kg   NaCl  as  an  independent  check  on  the  value  of 
pHCS).   The  solution  with  the  highest  concentration  of  NaCl 
CO. 06m)  corresponds  to  an  ionic  strength  of  0,16  mol  kg   . 
Figure  3  is  a  plot  of  pC^Tq^    from  equation  3.4  as  a  func- 
tion of  molality  of  NaCl  added  at  25  and  37°C.   The  inter- 
cepts, PtanYci^°  a"t  mci  =  °'  found  by  linear  regression  anal- 
ysis, were  7.520  and  7.499  at  25  and  37°C,  respectively, 
with  a  mean  deviation  slightly  greater  than  0.001.   Intro- 
duction of  the  pH  convention  to  equation  3.5  gave  values  of 
pHCS)  of  7.411  and  7,387  at  the  two  temperatures;  these  com- 
pare well  with  the  NBS  assigned  values.   Table  5  gives  the 
results  of  these  and  similar  measurements  involving  the  TES, 
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HEPES,  and  Tris  buffers  and  the  ph.osph.ate  buffer,  where  its 
contribution  to  the  total  ionic  strength  was  varied  by  ad- 
justing the  phosphate-Nad  ratio.   All  solutions  contained 
enough  Nad  to  bring  the  ionic  strength  to  0.16  mol  kg"1. 
Assuming  the  pH  convention  is  valid  at  this  ionic  strength 
(.above  its  intended  limit  of  0.1m),  pH(S>  values  were  as- 
signed to  these  buffer  solutions  based  on  equation  3,4.. 
They  are  listed  in  Table  5  with  the  corresponding  EMF's  and 
t  emper atur es . 

Measurements  in  a  cell  with  liquid  junction  have  been 
carried  out  on  the  same  solutions,  and  these  results  are 
listed  in  Table  6.   The  pH  values  listed  here  were  calculated 
from  the  EMF  values  and  equation  3.7,  referenced  to  the  EMF 
of  the  dilute  phosphate  buffer  and  its  assigned  pHCSl  value. 
The  values  of  pHCXl  calculated  in  this  manner  differ  appre- 
ciably from  those  assigned  by  the  conventional  approach. 
Both  sets  of  values  are  summarized  for  the  TES,  HEPES,  and 
Tris  buffers  in  the  second  and  last  columns  of  Table  7,   It 
is  discrepancies  of  this  sort  that  have  led  others29'69  to 
postulate  that  a  residual  liquid-junction  potential  exists 
between  dilute  and  higher  ionic  strength  buffers  when  meas- 
ured in  cells  with  liquid  junction. 


Figure  3.   Plot  of  pCa^)  (or  the  NBS  phosphate  "Mood 

buffer"  as  a  function  of  molality  of  added  NaCl 
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Isopiestic  Measurements 

NaCl-SrCl2  Mixtures 

Of  the  aqueous  systems  10M0C2  ,  where  M  =   Na,  K,  and 
N  =  Mg,  Ca,  Sr,  Ba,  only  the  system  NaCl-SrCl2  has  not  been 
studied  by  the  isopiestic  technique.   Thermodynamic  proper- 
ties of  this  system  are  of  interest  to  the  chemical  ocean- 
ographer,  since  both  components  are  present  in  seawuter. 

The  literature  values4  for  the  osmotic  coefficients  of 
strontium  chloride  in  water  are  based  upon  the  results  of 
two  isopiestic  studies,  one  using  KC170  and  the  other  CaCl  71 
as  isopiestic  reference.   Downes72  has  drawn  attention  to  the 
poor  agreement  between  values  derived  from  these  studies  and 
has  made  measurements  using  KC1  as  a  reference.   Measurements 
have  been  made  in  this  laboratory73  using  NaCl  as  a  reference 
up  to  about  3  mol  kg"1;  these  are  in  excellent  agreement  with 
Downes >  results.   Comparison  of  the  literature  values  with 
those  from  the  present  study  are  given  in  Table  8. 

Data  for  the  single  electrolytes,  NaCl  and  SrCl2 ,  have 
been  fitted  to  equations  4.6  and  4,8,  and  the  best-fit  param- 
eters have  been  obtained.74  These  parameters,  listed  in 
Table  9,  were  used  in  fitting  the  data  on  the  mixtures.   The 
compositions  of  the  mixtures  in  isopiestic  equilibrium  are 
given  in  Table  10,  along  with  the  molalities  of  the  NaCl 
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reference  solutions,  and  the  osmotic  coefficients  calculated 
from  equation  4,3.  Osmotic  coefficients  for  the  NaCl  refer- 
ence solutions  were  obtained  from  the  tables  of  Robinson  and 

4 
Stokes . 

The  osmotic  coefficients  of  the  mixtures  were  fitted  to 
the  Scatchard  neutral -'electrolyte  treatment  (equation  4.11); 
the  best-fit  parameters  and  standard  deviations  cr  (.$ )  are 
given  in  Table  11.   The  computer  program  allows  the  inclusion 
of  any  or  all  parameters  at  the  choice  of  the  operator.   This 
enables  one  to  see  the  significance  each  parameter  has  on 
the  standard  deviation  of  fit. 

The  same  procedure  has  been  applied  to  the  Pitzer  treat- 
ment (equation  4.12),  and  the  results  are  given  in  Table  11. 
It  can  readily  be  seen  that  the  Pitzer  treatment  for  the 
single  electrolytes  is  far  superior  to  that  of  Scatchard, 
in  that  these  equations  alone  fit  the  data  for  the  mixtures 
with  a  standard  deviation  a  (,<J> )  =  0.002  without  the  inclusion 
of  any  mixing  parameters,  while  the  Scatchard  treatment  re- 
quires at  least  one  mixing  parameter. 

NaC.l-Na2C0   Mixtures 

The  carbonate  system  in  seawater  is  important  from  the 
standpoint  of  acid-base  processes  which  influence  the  pH  and 
buffering  capacity  of  this  medium  as  well  as  the  solubility 
of  marine  carbonates.   It  is  of  interest  to  investigate  the 
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mixtures  NaCl-Na^  at  various  compositions  and  ionic 
strengths  to  elucidate  osmotic  coefficient  behavior  and 
obtain  thermodynamic  information  based  on  interactions  of 
these  electrolytes. 

Data  for  the  single  electrolytes  have  been  obtained,75 
and  the  useful  parameters  for  treatment  of  the  mixtures  are 
given  in  Table  12  along  with  the  standard  deviations  of  fit. 
The  isopiestic  molalities  and  osmotic  coefficients  of  the 
mixtures  are  listed  in  Table  13.   Blanks  in  the  table  are 
where  data  were  not  included  due  to  poor  precision  between 
duplicate  cups.   Results  of  fitting  the  mixtures  to  the 
Scatchard  and  Pitzer  treatments  are  given  in  Table  14. 
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Table  8.   Osmotic  Coefficients  in  Aqueous  SrCl„  Solutions 
at  25°C;  Comparison  of  Results 


SsrCl. 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

2.5 

3.0 

3.5 

4.0 

4,2 
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4>    CKCl   ref .  )    w      j>    (CaCl2    ref .  ) 


71 


0.8494 
0,8491 
0.8590 
0,8734 
0.8906 
0.9097 
0,9303 
0.9520 
0.9746 
0.9978 
1.0460 
1.0953 
1.1460 
1.1982 
1.2530 


0,8490 

0.8514 

0,8636 

0.8797 

0.8979 

0,9177 

0,9388 

0.9611 

0,9843 

1.0086 

1 , 0598 

1.1144 

1.1718 

1.2320 

1,2944 

1.4585 

1,6294 

1 . 80Q9 

1,9669 

2.0302 


<r 

<f>  CNaCl  ref.  1 

0.85Q 

0,8487 

0.854 

0,8501 

0,864 

0.8615 

0,880 

0,8768 

0 .  899 

0,8946 

0,918 

Q,ai40 

0.937 

0,9349 

0.959 

0,9570. 

0,983 

Q,98Q2 

1.009 

1,Q045 

1 ,  061 

1,Q558 

1,116 

1,1104 

1.173 

1,1674 

1.232 

1.2281 

1,292 

1.2905 

1.454 

lv4551 

1.631 

1.630Q 

1.802 

1 , 8010 

1,966 

1.9669 

2.0302 

Present  work. 
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Table  9,   Parameters-  for  Eqs.  4.6  and  4,8  Used  in  the 
Treatment  of  NaCl-SrCl2-H20  Mixtures 


Parameters  for  Eq. 

4.6 

S  -   -1,17202 

NaCla 

SrCl2b 

A 

1,4635 

1,60.11 

B 

0.041340 

0,031312 

C 

0.02083Q 

0.0095259 

D 

-0.0.016130 

-0.0Q05196Q 

E 

0.000043460 

Q,000Q13640 

I (max. ) 

6,0 

9.0 

a0j>) 

0 , O0Q4 

0,0012 

Parameters   for 

A®   =    -0. 

Eq. 
392 

4,8 

gC0) 

NaClC 
0.Q7669 

SrCl2b 
0.28994 

gCi) 

0,26461 

1.5795 

c* 

0.0012193 

-0.003755 

I  (max .  ) 

6.0 

6,0 

*(.*) 

0.0007 

0.0017 

Parameters  taken  from  ref.  74. 

Combined  fit  to  results  of  ref.  72  and  present  work 
CTable  8) . 


Parameters  obtained  by  fitting  smoothed  values  of  <j>„  „, 
from  ref.  4.  JNaL1 
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TableM1J^A?eSt:Fit  MixinS  Parameters  for  the  System 


NaCl(A)  -  SrCl2(B)  -  HgO 


*01 


0.0096 
0.0212 
0.0235 
0.0243 
0.0178 


Scatchard  neutral-electrolyte  treatment 


^02 


-0.0027 
-0.0040, 
-0.0038 
0.0002. 


^03 


0.0001, 


-0.0005 


*12 


-0.0014 
-0.0003, 


*13 


-0.0002, 


0.0050 
0.0017 
0.0011 
0.0011 

0.0005 
0.0005 


4 


Pitzer  treatment 


a,Sr 


0.0036 
-0.0096 
-0.0084 
-0.0076 


-Na,Sr 


0.0031 


0.0068 


-Na,Sr,Cl 


0.0036 
-0.0052 


0.0021 
0.0017 
0.0014 

0.0015 
0.0013 
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Table   12.      Parameters   for  Eq.    4.6   and  4.8  Used   in   the 
Treatment   of  NaCl-Na2C03-H20  Mixtures 

Parameters   for  Eq.    4.6 

S  =  -1.17202 

a 


A 
3 
C 
D 
E 
I(max) 
ffC4>) 


NaCl' 

1.4635 
0.041340 
0.020830 
-0.0016130 
0.000043460 
6.0 
0.0004 


Na2C03' 


1.4226 
-0.065322 

0.0084332 
-0.00028700 

0.0000018700 

9.0 

0.0020 


Parameters  for  Eq.  4.8 
A*  =  -0.392 


NaCl" 

NaoC0o 

— 2 — 3 

p(0) 

0.076690 

0.040822 

6(1) 

0.026461 

1.4679 

c* 

0.0012193 

0.004237- 

I (max) 

6.0 

9.0 

CJ($) 

0.0007 

0.0020 

Parameters  taken  from  ref.  74. 

Parameters  obtained  by  fitting  smoothed  values  of 
from  ref,  75. 


'Na2C03 


Parameters  obtained  by  fitting  smoothed  values  of  ± 
from  ref.  4.  ^NaCl 


Parameters  taken  from  ref.  75, 
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Table  14,   Best -Fit  Mixing  Parameters  for  the  System 
NaCl(A)  -  Na2C03CB)  -  H90 


Sea 


-0.1099 

-0.0956r 

-0.1744 

-0.1654 

-0.136. 


Scatchard  neutral-electrolyte  treatment 


*02 


-0.0036. 


0.0424. 


0.0197. 


-0.0243, 


^03 


-0.0061, 


0.00392. 


*12 


0,0180. 
0.0147 


—13 

c(i>) 

— 

0.047 

- 

0.0081 

- 

0,0081 

- 

0.0073 

_ 

0.0047 

3 


0,001723  0.0047 


-Na , Sr 

-0.0495, 
-0.0898- 
-0,1605 
-0.1607 


Pitzer  treatment 


-Na , Sr 


0.0103. 


-0.00254 


il» 


Na,Sr,Cl 


0.0503 
0.0548, 


0 


0.030 

0.012 

0.010 

0,0047 

0.0047 


CHAPTER  7 
DISCUSSION 


Electromotive  Force  Measurements 


Activity  Coefficient  of  HC1  in  Seawater 

For  the  HCl-MgCl2  system,  two  values  of  the  Harned 
slope,  a,    at  I  =  0.1  mol  kg-   and  25° C  have  been  reported 

previously.   That  of  Downes,    based  on  measurements  made 

77  7R 

by  Smith,    is  0.0521,  while  Khoo  et  al .    found  a  =  0.0326, 

This  is  a  rather  large  discrepancy.   The  value  obtained 
from  this  work,  a  =  0.0493,  tends  to  support  the  value 
reported  by  Downes. 

The  primary  purpose  of  this  study  was  to  derive  values 
of  a  at  ionic  strengths  corresponding  to  those- of  seawater 
of  salinities  20,  35,  and  45  /00,  and  thus  to  elucidate 
the  behavior  of  activity  coefficients  in  solutions  resem- 
bling seawater.   It  is  also  of  interest  to  evaluate  the 
trace  activity  coefficient  of  HC1  in  the  mixtures  and  com- 
pare it  to  earlier  measurements.   This  has  been  done  in 
two  ways.   First,  values  of  a  obtained  from  the  present 
work  involving  HCl-MgCl2  (Table   3)  were  combined  with 
those  of  an  earlier  study  of  HCl-NaCl,  weighting  the  two 
in  accord  with  the  relative  contributions  of  the  two  salts 
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to  the  total  ionic  strength  of  a  seawater  solution  con- 
taining only  those  two  salts  (equation  2.40).   In  another 

tr 
approach,  YHC1  was  calculated  from  the  Harned  slope  ob- 
tained from  HCl-NaCl-MgCl2  mixtures  at  I  =  0,6729  (Table 
3),   treating  NaCl  +  MgCl2  as  one  component  (equation  2.38), 
where  y  is  the  ionic  strength  contribution  of  the  two  salts 

to  the  total  ionic  strength.   In  Table  15,  the  values  of 

tr 
YHC1  stained  in  these  two  ways  are  compared  with  the 

values  obtained  by  Khoo  et  al .  from  EMF  measurements  of 
0.01m  HC1  in  synthetic  seawater.   Although  sulfate  was 
omitted,  the  salt  mixture  used  in  the  latter  study  approached 
closely  the  composition  of  natural  seawater,  being  composed 
of  NaCl,  MgCl2,  CaCl2,  and  KC1  in  the  correct  proportions. 
An  examination  of  the  table  shows  that  yH,  derived  from 
the  Harned  rule,  using  a  measured  for  HC1  in  the  salt  mix- 
ture composed  of  both  NaCl  and  MgCl2 ,  agrees  well  with  that 
calculated  by  equation  2,40  from  the  individual  values  of 

a  from  the  HCl-NaCl  and  HCl-MgCl2  systems.   Furthermore, 

tr 
it  is  evident  that  YHC1  is  nearly  the  same  in  a  mixture  of 

the  two  salts  as  in  a  mixture  of  the  same  ionic  strength 

but  containing  the  four  salts  predominant  in  seawater,   The 

average  difference  at  the  nine  temperatures  is  less  than 

0.001;  for  a  solution  of  NaCl  alone  (I  =  0.6729  mol  kg"1), 

tr 
^HCl  differs  from  that  in  seawater  on  the  average  by  about 

0,006,   Thus,  the  drop  in  YHC1  from  pure  HC1  to  seawater  is 
accounted  for  by  the  effect  of  the  two  major  salts,  NaCl 
and  MgCl2, 
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Table  15.  Comparison  of  the  Trace  Activity  Coefficient  of 
HC1  (Y^r)  in  NaCl-MgCl2  Mixtures  With  That  Determined 
Experimentally  in  Synthetic  Seawater  Csw) 


1  =  0. 

3809 

I  = 

0.6729 

1  =  0. 

8720 

CS  =  2C 

)°/00) 

CS  = 

35°/00) 

CS  =  4£ 

°/00). 

t,  °c 

Eq.  2.4Qa  swb 

Eq 

.  2,40 

Eq,  2.38 

c 
sw 

Eq.  2.40   sw 

5 

0.740 

0.740 

0 

741 

0.741 

0.741 

0.752 

0.752 

10 

0.738 

0.738 

0 

.739 

0.739 

0.738 

0,747 

0.748 

15 

0.739 

0.739 

0 

736 

0.736 

0.735 

0.747 

0.745 

20 

0.733 

0.737 

0 

734 

0,733 

0.732 

0.744 

0.741 

25 

0.730 

0.733 

0 

730 

0.730 

0.729 

0.739 

0,738 

30 

0.727 

0.730 

0 

726 

0.727 

0.725 

0.737 

0.734 

35 

0.725 

0.728 

0, 

722 

0.723 

0.722 

0,731 

0.730 

40 

0.722 

0.726 

0. 

718 

0.719 

0.718 

0.726 

0.726 

45 

0.718 

0.722 

0, 

713 

0.714 

0.714 

0.722 

0.722 

Calculated  by  Eq,  2.4Q  from  a  values  for  HCl-NaCl  mixtures 
(.reference  24).  and  for  HCl-MgCl2  mixtures  (present  work). 

b  _ 

From  experimental  data,  reference  25, 


c 


Calculated  by  Eq.  2.38  from  a  for  HCl-NaCl-MgCl-  mixtures 
yNaCl/MgCl0  =  1. 
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In  this  respect,  it  is  of  interest  to  consider  the 
results  in  light  of  Pitzer's  recent  treatment  concerning 
electrolyte  mixtures.   The  data  for  HC1  in  mixtures  of  NaCl 
and  MgCl2  have  been  fitted  to  the  Pitzer  equation  Cequation 
4.12)  in  an  attempt  to  estimate  the  pertinent  mixing  terms. 
The  parameters  8,^)  and  6?^)  for  the  single  electrolytes 
were  taken  from  Pitzer's  tabulated  values.48  Strictly, 
the  CMX  terms  should  have  been  included.   However,  they 
were  unavailable  and  were  taken  to  be  zero.   The  best -fit 
values  of  interaction  terms  9^,  ^,  and  ^  were  derived 
by  curve-fitting.   The  results  are  summarized  in  Table  16. 
It  is  evident  from  the  standard  deviations  given  in  the 
last  column  that  the  properties  of  the  mixtures  are  not 
well  accommodated  by  8(0)  and  S(1)  for  the  individual 
electrolytes  alone.   Inclusion  of  the  cation-cation  mixing 
term,  G,  produces  considerable  improvement.   The  parameter, 
Gf  ,    which  represents  a  dependence  of  0  on  ionic  strength, 
appears  to  be  of  little  importance  in  fitting  the  data. 
The  ternary  interaction  term,  i|> ,  produces  some  improvement; 
however,  its  inclusion  would  result  in  a  term  quadratic  in 
yMgCl   in  Pi"tzer'sr  equation  for  the  activity  coefficients 
in  the  mixture,   Since  the  data  obey  Harned's  rule  without 
such  a  term,  the  necessity  of  0  is  questionable. 

The  Harned  slope  appears  to  be  linearly  dependent  on 
temperature,  as  shown  by  the  slope  da/dT  given  at  the  foot 
of  Table  15.   Since  -ayBI  =  log  (Yhc1/y£c1),  the  quantity 
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Table  16        Mixing  Parameters    (Pitzer  Treatment)    for  the 
System  HCl-MgCl2-H20 


§H,Mg 

-H,Mg 

^H,Mg,Cl 

103s 

a 

- 

- 

- 

30.9 

-0.2904 

- 

- 

5.8 

-0.3430 

0,2728 

- 

4.1 

-0.4357 

- 

0,3843 

2.9 

-0.4651 

-0.1401 

0.5337 

2.8 

Standard  deviation  of  log  y 


81 
da/dT  may  be  related  to  the  difference  in  partial  molal 
heat  content  of  HC1  in  its  pure  solution  and  in  a  mixture 
with  MgCl2  by  the  equation: 

-4.606RT2yBI(da/dT)  =  2RT2d_  In  (Yhc1/yJc1) 

=  LHCl^MgCl2  '  °)  "  W^MgCl^ 

(7.1) 
Thus,  the  partial  molal  enthalpy  appears  to  be  linearly 
dependent  on  yB>  and  in  this  respect  an  analog  of  Harned's 
rule  exists,  at  least  for  5°C  <  T  <  45°C  and  0.1  <  I  <  0.9. 

pH  Standardization  in  Seawater 

Comparison  of  results  (Table  4)  for  prn^X)  obtained  by 
the  operational  definition  with  pmH(S)  values  calculated  from 
equation  3.15  shows  that  a  discrepancy  does  exist  when  salin- 
ity is  varied.   At  a  salinity  of  30°/00,  for  example,  this 
amounts  to  an  average  of  0.007  pH  unit  higher  for  the  prn^X) 
values,  while  at  S  =  40°/00,  the  operational  values  are,  on 
the  average,  0.013  pH  unit  lower  than  the  corresponding 
pmH(S)  values.   Within  a  given  salinity,  however,  pmH(X) 
values  are  essentially  identical  regardless  of  the  Tris 
buffer  concentration,  indicating  a  stabilization  of  the 
residual  liquid-junction  potential  by  the  seawater  medium. 

It  is  interesting  to  note  that  when  the  operational 
definition  is  employed  using  Ex  with  Eg  and  pH(S)  of  the  NBS 
1:3.5  phosphate  buffer,  the  resulting  pH^  values  are  very 
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close  (0.007  pH  unit  on  the  average)  to  the  corresponding 

pm^CS)  values,  even  though  there  is  a  large  difference  in 
n 

ionic  strength  between  the  seawater  solutions  and  the  dilute 
(I  =  0.1m)  phosphate  buffer. 

Since  pH(X)  differs  from  pH(S)  (=  -log  nigYg)  by  the 
liquid- junction  term,  E./k,  the  relationship  between  pHNBg 
and  pmf.(S)  can  be  written  as: 

pmgCS)  =  pHNBg  +  log  yH  +  Ej/k  (7.2) 

The  close  agreement  between  pHNBg  and  pmH(S)  must  then  be 
attributed  to  a  cancellation  of  the  liquid-junction  poten- 


tial and  activity  coefficient  terms.   Robinson  and  Bates 
have  recently  estimated  yc1  =  0.639  in  seawater  of  35  /0C 
salinity,  based  on  an  hydration  number  method  developed 
earlier.67  Using  yH-.=  0.729  at  25°  C,  one  obtains  0.832 
Trr  (=  Ywi/Ym)'   Using  this  value  in  equation  7.2  with 
pnu.CS)  -  pH„DC,  =  -0.007  at  35  /00  salinity,  the  residual 
liquid- junction  potential  is  estimated  to  be  about  4.3mV. 


pH  Standardization  in  Clinical  Media 

As  pointed  out  in  Chapter  6,  discrepancies  between 
column  2  and  the  last  column  of  Table  7  are  thought  to  arise 
as  the  result  of  an  appreciable  residual  liquid-junction 
potential  caused  by  the  difference  in  ionic  strength  between 
the  relatively  dilute  NBS  phosphate  "blood  buffer"  and  the 
higher  ionic  strength  (I  =  0.16m)  Tris,  HEPES,  and  TES 
buffers. 
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The  pH(S)  values  obtained  for  the  phosphate  buffer  + 
NaCl  by  measurements  in  a  cell  without  liquid  junction  in  all 
likelihood  lie  close  to  the  NBS  standard  scale,  based  on  the 
same  pH  convention.   The  reasoning  is  as  follows. 

Calculations  from  equation  3.3  for  the  logarithm  of  the 
activity  coefficient  of  chloride  ion  at  25° C  give  -0.110  at 
I  =  0.01m  and  -0.128  at  I  =  0.16m.   These  values  are  nearly 
the  same  as  the  logarithm  of  the  mean  activity  coefficient 
(log  y±)  of  NaCl  at  I  =  0.1m  (-0.109)  and  at  I  =  0.16m 
(-0.126).   At  37° C,  log  ycl  is  -0.130  at  I  =  0.16m,  while 
log  y+  is  -0.128.   Furthermore,  the  hydration  convention79 
yields  -0.130  for  log  yc1  in  a  solution  of  NaCl  of  molality 
0.16  mol  kg"1. 

One  expects  that  log  y±  and  log  yc1  will  vary  linearly 
as  phosphate  is  added  to  a  solution  of  NaCl  maintaining  a 
constant  ionic  strength.0  To  estimate  the  magnitude  of 
this  change,  the  sodium  glass  electrode  was  employed  in  a 
cell  of  the  type: 

NalSE | Phosphate  buffer  +  NaCl|AgCl;  Ag        (7.3) 
and  the  following  EMF  data  were  obtained: 
ionic  strength  E(V) 

phosphate     NaCl      25° C       37° C 

0  0.16      0.0432      0.0465 

0-04        0.12      0.0540      0.0574 
0-10        0.06      0.0767      0.0814 
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The  standard  potential  of  cell  7.3  was  derived  knowing 
log  y+  in  the  0.16m  NaCl  solution,  and  the  activity  coeffi- 
cients in  the  mixtures  were  calculated  by: 

-  log  Y±  =  (E  "  -  )  +  i    log  (mNamcl)       (7.4) 

The  results,  plotted  in  Figure  4,  are  not  of  high  accuracy 
but  do  demonstrate  a  linear  increase  of -log  y+  when  NaCl  is 
replaced  by  buffer.   This  increase  amounts  to  only  0.001  in 
log  y+  for  each  increment  of  0.01  in  the  ionic  strength  con- 
tribution of  phosphate.   Therefore  any  discrepancy  in  calcu- 
lation of  -log  yc1  by  the  convention  for  a  solution  of 
phosphate  (I  =  0.06)  and  NaCl  (I  =  0.1)  would  not  amount  to 
more  than  0.006,  small  when  compared  with  the  difference  of 
0.042  at  25° C  and  0.030  at  37° C  in  pH  values  derived  from 
cells  with  and  without  liquid  junction  (see  Tables  5  and  6). 

Having  values  of  pH(S)  for  buffer  solutions  with  NaCl 
(last  column,  Table  7)  that  appear  to  lie  close  to  the  NBS 
standard  scale,  it  was  possible  to  examine  the  effect  of 
ionic  strength  and  composition  of  the  phosphate  standard  on 
the  pH(X)  values  determined  for  the  TES,  HEPES,  and  Tris 
buffers  from  the  cell  with  liquid  junction.   By  taking 
values  of  E  and  pH(S)  for  the  reference  phosphate-NaCl 
solutions  from  Tables  6  and  5  respectively  and  the  values  of 
Ex  for  the  other  buffers  (Table  6)  and  applying  equation  3.7, 
the  pH(X)  values  listed  in  Table  7  were  calculated;  each 
column  heading  being  the  reference  buffer  composition.   The 


Figure  4.   Variation  of  the  logarithm  ^  *fc 

sltion  at Is'and  3^?'  aS  '  fUnCtlon  of  <"»<Po- 
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results  are  shown  graphically  in  Figure  5.   The  horizontal 
dotted  lines  indicate  the  values  of  pH(S)  assigned  to  the 
respective  buffers. 

It  is  clear  from  Figure  5  that  the  greater  part  of  the 
discrepancy  between  the  operational  pH(X)  and  the  conventional 
activity  pH(S)  is  removed  simply  by  matching  the  ionic 
strength  of  the  standard  to  that  of  the  unknown  sample.   By 
increasing  the  NaCl  molality  further  at  the  expense  of  phos- 
phate such  that  I  (phosphate)  =  0.06m  and  I  (NaCl)  =  0.10m..  the 
residual  liquid-junction  potential  is  completely  eliminated 
for  the  TES  and  HEPES  buffers  although  a  difference  of  about 
0.9mV  (0.015  pH  unit)  persists  for  the  Tris  buffer. 

Estimates  of  liquid-junction  potentials  at  the  junction 

buffer  |  |3.5M  KC1  have  been  made  by  use  of  the  Henderson 

equation.0   Using  available  limiting  ionic  conductances  at 

25° C  and  estimating  those  of  HPO  2~,  TES-,  HEPES-,  and 

TrisH  ,  the  following  values  have  been  obtained: 

buffer  I,  m    Ej ,  mV 

1:3.5  phosphate 
P(0.10)  +  NaCl(0.06) 
P(0.08)  +  NaCl(0.08) 
P(0.06)  +  NaCl (0.10) 
P(0.04)  +  NaCl (0.12) 
1:1  TES (0.04)  +  NaCl (0.12) 
1:2  TES(0.04)  +  NaCl (0.12) 
1:2  HEPES(0.04)  +  NaCl(0.12) 
1:3  Tris (0.05)  +  NaCl (0.11) 

It  is  evident  that  the  residual  liquid-junction  potential 

when  the  blood  buffer  (I  =0.1m)  is  replaced  by  a  TES  buffer 

(I  =  0.16m),  for  example,  should  not  exceed  0.6mV  (0.01  pH 

unit),  whereas  the  observed  differences  are  three  times  this 
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Figure  5.   Effect  of  the  ionic  strength  and  composition  of 
the  reference  standard  on  the  operational  pH 
values  of  TES,  HEPES,  and  Tris  buffer  solutions 
of  ionic  strength  0.16m  at  37°  C 


89 


7.67- 


X 

a. 


720 


< > 

■  0.04  P/QJ2  NaCI 


•  0.06     /O.J0 

A  0.08     /0.08 

0.10     /0.06 


■J— '{without  I.].) 


a«  QJ6 

Ionic  Strength 


90 

figure.   However,  the  calculations  do  lead  to  the  conclusion 

that  the  potential  difference  should  be  small  when  a  buffer 

0,06m  in  phosphate  and  0.1m  in  NaCl  is  measured  against  the 

TES  and  HEPES  buffers  of  the  same  ionic  strength. 

The  low  value  of  E.  calculated  for  the  Tris  +  NaCl 
3 

buffers  compared  with  the  other  buffers  of  the  same  ionic 
strength  is  the  result  of  a  larger  value  of  the  difference 
in  the  limiting  ionic  conductances,  Zc  \°    -  Zc  X°      which,  in 
turn,  is  caused  by  the  fact  that  the  mobility  of  the  TrisH 
cation  is  lower  than  those  of  Na   and  K  ,  the  cations  of 
which  the  other  buffers  are  composed.   The  three  hydroxyl 
groups  on  the  Tris  structure  and  their  capacity  for  hydrogen- 
bonding  to  water  support  this  suggestion. 

Isopiestic  Measurements 

Scatchard  Treatment 

In  the  simplest  case,  the  osmotic  coefficient  of  a 
solution  composed  of  1:1  and  a  1:2  or  2:1  electrolyte  with 
a  common  ion  will  be  a  simple  function  of  the  contributions 
of  single  electrolyte  properties  as  described  by  the  trunca- 
ted form  of  equation  4.11: 

(1  +  yA}*MIX  =  2*AyA  +  *ByB  (?'5) 

This  type  of  behavior  would  be  observed  if  only  binary  inter- 
actions between  ions  of  opposite  charge  were  present,  as  first 

13 
postulated  by  Guggenheim.     In  reality,  however,  particularly 
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at  higher  ionic  strengths,  it  becomes  necessary  to  consider 
binary  interactions  between  like  charged  ions  as  well  as 
additional  higher  order  (e.g.,  ternary)  interactions.   This 
forces  the  inclusion  of  additional  terms  to  equation  7.5  if 
one  wishes  to  describe  the  data  mathematically.   Scatchard's 
general  equation  for  the  systems  studied  here  may  be  written 
as : 

(1  +  yA^MIX  =  2*AyA  +  *ByB  +  yAyBe0'  +  yAyB(yA  "  yB)6l  ^ '^ 
where  6Q  =  b^I  +  b^I2  +  b^I3  +  .  .  .  and  61  =  b^I2  +  b^I3 
+  . . .   The  functional  dependence  on  m^  and  mB  of  the  terms 
in  which  the  mixing  parameters,  b .  . ,  reside  has  led  Robinson80 
to  identify  each  parameter  with  a  particular  type  of  ionic 
interaction,  none  of  which  is  characteristic  of  only  the 
single  salt  solutions.   The  bQ1  term  is  identified  with 
binary  interactions  between  ions  of  like  charge,  b  9  with 
ternary  interactions  involving  two  like-charged  ions  and 
one  of  opposite  charge,  and  b12  with  ternary  interactions 
among  ions  all  of  like-charge.   The  terms  bn„  and  b10  are 
sometimes  associated  with  quaternary  interactions  but  are 
most  likely  only  "data-fitting"  parameters.   The  inclusion 

of  any  additional  terms  is  superfluous. .  The  Oak  Ridge  com- 

81 
puter  program  used  in  this  work   allows  the  evaluation  of 

the  first  five  parameters. 

Examination  of  the  top  of  Tables  11  and  14  shows  that 

while  the  inclusion  of  all  five  of  the  above  parameters 

gives  the  best  fit  to  the  experimental  results,  the  same 
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agreement  can  be  obtained  using  only  b01,  b„2,  and  b10. 
The  relative  magnitude  oi  the  parameters  between  the  NaCl- 
SrCl2  and  the  NaCl-Na2C03  studies,  as  well  as  the  inferior 
fit  of  the  latter  data,  indicate  much  more  complex  behavior 
in  the  system  involving  carbonate.   This  is  reinforced  by  a 

comparison  of  Figures  6  and  7,  in  which  a  plot  of  d>,,T-  as  a 

MIX 

function  of  I  has  been  made  for  the  two  systems.   It  can 
readily  be  seen  that  due  to  the  large  difference  in  <£°  for 
the  pure  solutions  NaCl  and  Na2C03,  $    is  highly  dependent 
on  yNa  CO  at  h^h   ionic  strength  while  this  is  not  the  case 
in  the  NaCl-SrCl2  system. 

The  free  energy  of  a  solution  of  a  single  electrolyte 
is  made  up  of  the  free  energy  of  the  components  in  their 
standard  state  and  the  ideal  free  energy  which  is  a  function 
of  the  concentration  of  solute.   Any  additional  contribution 
is  referred  to  as  excess  free  energy  and  is  made  up  of  Debye- 
Hiickel  forces  and  any  additional  interaction  terms.   It  is 
this  excess  free  energy  which  is  estimated  by  the  activity 
and  osmotic  coefficients. 

When  two  solutions  are  mixed,  there  arises  an  excess 
free  energy  of  mixing  above  and  beyond  that  present  in  each 
respective  solution,  and  it  is  this  quantity  for  which  the 
interaction  parameters  may  have  some  significance.   It  is  of 
interest,  therefore,  to  calculate  this  excess  free  energy  of 
mixing  for  our  two  systems.   From  the  last  two  terms  on  the 


Figure  6.   Plot  of  ^  as  a  function  of  ionic  strength  in 
the  system  NaCl(A)  -  SrCl2(B)  -  Ho0  for 
ionic  strength  fractions,  y 
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Figure  7.   piot  of  ^  as  a  ^.^  ^    .^  ^^ 

«  the  system  NaCl(A)  -  Na  CO  (B)  -  HO  for 

various  ionic  strength  fractions   y 
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right  of  equation  7.6,  the  excess  free  energy  of  mixing 
per  kilogram  of  solvent  at  constant  ionic  strength  may  be 


73 
given  as : 


AG? 


MIX  =  RTIVB[B0  +  B1^A  "  V]  <7'7) 

where  BQ  =  bQ1I^+  (l/2bQ2)I2  +  (l/3b03)I3  and  B±   =    (l/2b12)I2 
+  (l/3b13)I3. 

t? 
The  AGMIX  for  each  system  has  been  calculated  for  the 

case  yA  =  yB  =  0.5  at  three  ionic  strengths  using  only  the 
bQ1  and  bQ2  values  listed  in  Tables  11  and  14.   The  results 
for  NaCl-SrCl2  are  listed  in  Table  17  with  other  members  of 
the  series  for  comparison.   The  NaCl-Na2C03  results  are  com- 
pared with  the  results  of  NaCl-Na2S04  from  another  study52 
(Table  17).   It  can  be  seen  that  there  is  a  trend  in  AG^, 
with  atomic  number  within  the  alkaline-earth  group,  and  that 
this  quantity  increases  with  ionic  strength  as  expected. 
From  the  results  of  the  NaCl-SrCl2  system,  it  would  appear 
that  the  relatively  slight  AGBJX  is  destabilizing  or  unfa- 
vorable due  to  the  positive  sign.   This  is  accounted  for 
solely  by  the  bQ1  term  and  implies  that  a  repulsive  energy 
arises  from  new  interactions  between  Na+  and  Sr++  which  are 

greater  than  the  sum  of  Cl~  -  Cl~,  Na+  -  Na+,  and  Sr+2  - 

+2 

Sr   repulsive  energies  that  no   doubt  exist  in  the  single 

electrolyte  solutions.   From  the  sign  of  the  bno  and  b10 
terms,  however,  it  appears  that  there  is  a  net  attractive 
energy  gained  among  ternary  interactions. 
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Table  17,   Excess  Free  Energy  of  Mixing  (J  kg"1) • 
Comparison  of  Results 


NaCl  -  MgCl2  -  H2055 


NaCl  -  CaCl0  -  H„0 
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NaCl  -  SrCl2  -  H20 


2    "2V 

81 


NaCl  -  BaCl2  -  HO 


NaCl  -  Na2C03 
NaCl  -  Na2S0452 


1 

3 

5 

35 

250 

578 

17 

130 

318 

12 

95 

222 

16 

46 

-11 

-96 

-758 

-1799 

-38 

-318 

-845 
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In  the  NaCl-NapCO„  system,  it  is  evident  from  the  nega- 
tive  AG„_^.  that  there  is  an  overall  stabilizing  effect  upon 
mixing.   The  sign  and  magnitude  of  b--  indicate  a  net 

attractive  energy  gained  in  new  interactions  between  Cl~ 

_2 

and  C0„  .   The  sign  of  bQ2  and  b-2  indicate  unfavorable 

ternary  interactions. 


Pitzer  Treatment 

In  the  lower  half  of  Tables  11  and  14  are  the  results 
of  the  Pitzer  treatment  on  the  two  systems  studied.   One 
recalls,  from  equation  4.12,  that  0  is  an  interaction  param- 
eter representing  additional  interactions  between  ions  of 
like  charge,  £'  is  its  possible  dependence  on  ionic  strength. 
and  i|>  corresponds  to  new  ternary  interactions  between  two 
ions  of  like  charge  and  one  of  opposite  charge.   It  can  be 
seen  from  the  results  of  the  NaCl-SrCl?  system  that  the 
Pitzer  treatment  fits  the  data  quite  well  using  only  infor- 
mation derived  from  the  single  electrolytes.   The  addition 
of  the  0  mixing  parameter  offers  some  improvement,  and  its 
sign  is  positive,  analogous  to  the  Scat chard  bn1  parameter. 
The  results  from  the  inclusion  of  additional  mixing  param- 
eters  are  ambiguous,  particularly  since  the  0  parameter 
changes  sign  and  magnitude  drastically  with  only  slight 
improvement  in  the  standard  deviation  of  fit.   Without  the 
need  for  parameters  other  than  0,  the  significance  of  the 
numbers  below  line  2  (lower  half,  Table  11),  beyond  being 
purely  data-fitting  coefficients,  is  questionable. 
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The  results  of  the  NaCl-Na2C03  system  indicate  that 
mixing  parameters  are  necessary  to  fit  the  data  reasonably, 
although  the  single  electrolyte  data  alone  produce  a  better 
fit  than  the  Scatchard  Treatment.   Inclusion  of  9  improves 
the  fit  substantially;  however,  it  is  necessary  to  include 
i>   to  attain  the  fit  given  by  the  Scatchard  equation  using 
b01'  b02 '  and  b12'   It:  is  obvious  that  there  is  no  signifi- 
cant dependence  of  9  on  ionic  strength,  hence  no  need  for 
9'.   The  signs  of  9  and  ty    are  consistent  with  Scatchard's 
bQ^  and  b  2  for  the  same  system. 

A  computer  program  has  been  written  for  the  calculation 
of  the  activity  coefficients  of  the  individual  solutes  from 
the  appropriate  form  of  the  Pitzer  equation,82  (Appendix) 
utilizing  the  parameters  from  both  the  single  electrolytes 
and  the  mixtures.   For  the  NaCl-SrCl2  system,  only  the  9 
mixing  term  was  used,  and  the  results  are  given  in  Figure  8. 
The  values  for  Y^aC1  were  taken  from  Robinson  and  Stokes4 
while  values  for  YjrC1   were  calculated  from  osmotic  coef- 
ficient data  on  the  single  electrolytes.   The  positive 
change  in  excess  free  energy  upon  mixing  is  easily  seen  as 
an  increase  in  the  mean  activity  coefficient  of  SrCl^  which 
far  outweighs  the  decrease  in  YNaC1  exhibited  at  high  ionic 
strength.   The  positive  sign  of  the  9  term  in  the  Pitzer 
equation  and  of  the  bQ1  term  in  the  Scatchard  equation  lead 
one  to  suspect  that  it  is  the  strontium  ion  which  is  so 
highly  affected  upon  mixing.   This  is  reasonable,  since  no 
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new  interactions  arise  between  anions  and  cations.   It  is 
interesting  that  the  sodium  ion  is  barely  affected,  and  only 
at  high  ionic  strengths.   Harned-type  plots  are  given  in 
Figure  9  for  this  system,  and  log  y+  is  seen  to  be  a  linear 
function  of  y  for  both  electrolytes  at  all  ionic  strengths. 
This  is  to  be  expected,  since  the  data  were  fitted--  adequately 
without  any  higher  order  terms  in  y,  and  0  alone  was  used  to 
generate  the  activity  coefficients. 

Figure  10  shows  the  activity  coefficients  in  the  mix- 
ture NaCl-Na2C0„.   The  negative  excess  free  energy  of  mixing 
is  evident  here  from  a  large  decrease  in  y„  „,  only  slightly 
compensated  for  by  the  increase  in  y^t  cq  .   The  sign  and 

_  9 

magnitude  of  the  mixing  parameters  point  to  the  CI   -  C0„ 
interaction  as  the  primary  cause  of  the  observed  behavior. 
It  can  readily  be  seen  also  that  addition  of  a  large  amount 
of  Na2C0„  to  a  smaller  amount  of  NaCl  will  have  a  greater 
effect  on  AG„,„  than  the  reverse  operation.   For  example, 

in  a  situation  where  y.  =  0.15  and  yR  =  0.85  at  I  =  3.0m, 

E  —1 

the  calculated  value  of  AG,.TV  is  -439  J  kg   ,  whereas  if 

MIX 

yA  =  0.15  and  yB  =  0.85,  AG^IX  =  -332  J  kg-1.   The  Harned 
plots  (Figure  11)  for  this  system  introduce  two  important 
points.   First,  all  the  plots  are  linear  except  for  log 
YN  ri  at  high  ionic  strength.   This  gives  evidence  that  it 
is  to  account  for  the  behavior  of  YNaC1  that  we  must  intro- 
duce the  higher  order  #  term.   Secondly,  it  appears  that  at 
lower  ionic  strengths,  those  which  apply  to  seawater,  the 
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data  might  adequately  be  fitted  with  the  0  term  alone.   In 
this  respect,  they  obey  Harned's  Rule.   The  intercept  of 
log  yb  at  yB  =  0  (yA  =  1)  at  an  ionic  strength  of  0.6m  gives 
for  the  trace  activity  coefficient  of  NaoC0„,  ytr  =  0.417. 
Using  the  estimate  of  0.709  for  yNa+  at  I  =  0.6m  made  by 
Robinson  and  Bates67  one  obtains  0.144  for  Yrn_2  which 
agrees  closely  with  that  estimated  by  Robinson.75 


Figure  8.   Mean  activity  coefficients  of  NaCl  and  SrCl 

as  a  function  of  ionic  strength  in  the  system 
NaCl (A)  -  SrCl2(B)  -H20  for  various  ionic 

strength  fractions,  y   and  y 

A       R 
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Figure  9.   Earned  plots  for  NaCl  and  SrCl2,  -log  Y+  vs 

-2(B)  -  H2( 


7B,  in  the  system  NaCl(A)  -  SrCl0(B)  -HO 
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Figure  10.   Mean  activity  coefficients  of  NaCl  and  Na„C0„ 
as  a  function  of  ionic  strength  in  the  system 
NaCl(A)  -  Na2C03(B)  -  H20  for  various  ionic 
strength  fractions,  y.  and  yR 
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Figure  11.   Harned  plots  for  NaCl  and  ^co   -log  Y  vs 
yB,  in  the  system  tfaCl(A)  -  Na2COg(B)  -HO 
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CHAPTER  8 
CONCLUSIONS 


The  primary  purpose  of  this  research  was  to  investigate 
the  factors  which  affect  the  thermodynamic  activity  of  indi- 
vidual species  in  complex  saline  mixtures.   By  utilizing  a 
variety  of  experimental  techniques  on  relatively  simple 
systems  of  which  complex  mixtures  (e.g.,  seawater)  are  com- 
posed, information  can  be  obtained  as  to  the  specific  influ- 
ences on  electrolyte  behavior. 

When  the  logarithm  of  the  mean  activity  coefficient  of 
HC1  is  plotted  as  a  function  of  yM_cl  ,  at  ionic  strengths 
corresponding  to  seawater,  a  straight  line  is  obtained. 
When  the  slope  of  this  line  ("Harned  Slope")  is  combined 
with  the  Harned  Slope  from  a  similar  plot  vs.  y„  _.,  weight- 
ing the  two  as  to  their  molal  ratio  in  natural  seawater, 
log  YHC1  beiiaves  iust  as  it  has  been  found  to  do  in  total 
seawater  excluding  sulfate.   In  light  of  recent  theoretical 
predictions  concerning  electrolyte  mixtures,  the  linear 
decrease  in  log  YHC1  appears  to  be  due  to  interactions 
among  the  cations,  H+  -  Na+  -  Mg+2 . 

An  adequate  knowledge  of  the  behavior  of  H  in  media 
of  this  sort  allows  the  formulation  of  a  pH  scale  for  the 
standardization  of  pH  measurements  in  seawater.   With  a 
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tr 
value  for  yhc1  and  a  reliable  estimate  of  yhc1  in  seawater, 

a  scale  of  pn^  in  addition  to  paH  is  feasible.   Further, 
examination  of  the  residual  liquid-junction  effects  which 
occur  when  a  comparison  is  made  of  measurements  in  cells 
with  liquid  junction  allows  the  estimation  of  the  error 
which  would  result  from  pH  measurement  in  seawater  of  a 
particular  salinity  when  the  glass  electrode  is  standard- 
ized with  a  more  dilute  NBS  phosphate  buffer.   This  residual 
liquid-junction  potential  has  been  estimated  to  be  about 
4.3mV  between  the  NBS  1:3.5  phosphate  buffer  (I  ■  0.1m)  and 
seawater  of  35°/00  salinity  (I  =  0.67)  at  25°C.   This  amounts 
to  about  0.07  pH  unit. 

It  has  been  shown  that  residual  liquid-junction  poten- 
tials, which  can  cause  errors  of  up  to  0.05  pH  unit  in  pH 
measurements  in  clinical  media,  can  be  eliminated  by  choice 
of  a  reference  buffer  where  phosphate  and  NaCl  contribute 
0.06  and  0.10m,  respectively,  to  the  total  ionic  strength 
(I  =  0.16m,  pH(S)  =  7.297  at  37°C).   This  pH(S)  value  lies 
close  to  the  NBS  primary  standard  scale,  and  thus  its  adop- 
tion as  a  primary  standard  is  practical.   Secondary  buffers 
of  TES  and  HEPES  with  pH  values  above  and  below  that  of 
blood  plasma  have  been  proposed. 

Finally,  isopiestic  studies  have  been  made  on  two  sys- 
tems involving  components  present  in  seawater,  in  an  attempt 
to  shed  light  on  the  activity  effects  brought  about  by 
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mixing.   Although  current  mixed  electrolyte  theories  allow 
one  to  pinpoint  the  ions  which  have  a  primary  effect  on  the 
free  energy  of  the  system,  much  work  remains  to  be  done  to 
elucidate  the  nature  of  these  interactions.   Speaking  only 
in  terms  of  attractions  and  repulsions  is  an  oversimplifi- 
cation at  best.   Trends  exhibited  in  electrolyte  behavior 
among  elements  in  a  periodic  family,  for  example,  lead  one 
to  suspect  that  ionic  charge  is  only  one  factor  and  that 
ionic  size  and/or  hydration  capabilities  are  responsible  for 
more  specific  phenomena.   Debye  and  Huckel  first  noticed  the 
need  for  an  ion-size  parameter,  and  this  provided  the  founda- 
tion for  all  current  treatments.   The  success  of  current 
theory  in  describing  electrolyte  behavior  is  formidable,  and 
continued  work  can  only  lead  to  more  useful  descriptions  of 
activity  effects  and,  ideally,  a  unified  theory  which  will 
provide  a  greater  understanding  of  the  nature  of  aqueous 
electrolyte  behavior  as  a  whole. 
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